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Die angenäherte Berechnung der Fourier-Koeffizienten 
einer periodischen Funktion 


Von ANDREI RIPIANU 


In dieser Arbeit wird versucht, eine Methode zu finden, mit der die Genauigkeit der Näherungswerte der 
Fourierkoeffizienten, die man bei der Unterteilung der Periode in s gleiche T’eile mittels eines der bekannten 
Rechenverfahren erhält, gesteigert werden kann. Zu diesem Zwecke suchen wir nach mathematischen Bezie- 
hungen zwischen den Näherungswerten der Fourierkoeffizienten im Falle der Unterteilung der Periode in s 
gleiche Teile und im Falle ihrer Unterteilung in ks gleiche Teile, ers k eine positive ganze Zahl ist. In der 
vorliegenden Arbeit wird der Fall k = p" behandelt. 


An attempt is made at finding a method that would allow to increase the accuracy of the approximate values 
of Fourier coefficients oblained by means of one of the well-known methods of calculation when the period is 
divided in s equal parts. To this end, mathematical relations are looked for that connect the approximate values 
of Fourier coefficients pertaining to a division of the period in s and ks equal parts, respectively, k being a 
positive integer. In the present paper, the case k = p" is discussed. 


B pa6oTe uımerca MeToA AA HAXO’KIEeHHA Ö0JIee TOYHBIX IIPHOJLFKÖHHEIX 3HAYeHNÜ Koahhu- 
IIMEeHTOB Dypbe, KOTOPbIe IIOAIY4YAIOTCH H3BECTHEIMH METONaMNM pasOneHuA MepnoNa Ha s PAaB- 
HbIx yacreh. ln 3Tof Mesim yCTaHABANBAITCH COOTHONIEHUA Meky HPHÖJNZRÖHHBIMN 3HA- 
YeHHUAMH KOAPPHIMEHTOB Dypbe B cıIyyae pasöneHnsn ITepmona Ha s PABHBIX YacTeüi U B caIyyae 


ee Ha k s PABHbIX yacTeü, k — HaTypaıbHoe yncno. B BEIONS PaccMmaTpuBaeTca cIy- 
yaük= 


Man betrachtet die periodische Funktion p(u). Man kann annehmen, daß die Periode gleich 
2 ist, daman dies immer durch einen passenden Wechsel der Veränderlichen verwirklichen kann. 


Die Entwicklung in eine FOoURIER-Reihe lautet 


o(u) = + I (a, cosvu + b,sinv u) 
ai 


mit 
2n 


27 
1 
a, =; [rw cosyudu, b, = — [w sinv udu, 1 U re 
IT IT 


1) 


Mit der angenäherten Quadraturformel: 


2 
1 22 an 
a ERW ER en 
- | Ku) du z 6a) [a 
0 
erhält man: 
(0) = Krtia)osvis, | u Ey eye a ah 
Sj-el \S/ S j=l 
Wir definieren: 
ee | 5, ® \ 
| a o(ia+Pß)cosvja, 8 pie EBsinvia sh Bl 
BER ß \ m (B) (8) 
Wir nehmen uns vor, | 5 und 4 als Funktion von | ) und | ) auszudrücken. 
Bi7S Dass Ss S 


p bezeichnet eine ganze positive Zahl, während ß eine Verschiebung des Integrationsbereiches 
(0, 2%) entlang der O u-Achse bezeichnet, die wir später festlegen werden (siehe Bild). 
32 


VER 


Q 
8 
8 


Man hat dann 


f ee n- 

und, weil Sa £ 

| - y(u) = RO De 

E ist, R i e 
Plprs +) al = 90 &,) (a = L, 2) Era". “ = f 


Aus (2) und (4) folgt: | ; 


(ß) m—175 
3 )- zm A & "HD im + B) 608» Piom 


+3 = Dit dm +Rlosr Bit, 


aus (3): 
y j P &m = &m—1» 
Daher ist 


») pm—1z 


a, (B) 1 
A Pr p a 5 & PU &m—ı + ß) os ym-ı 


el f j pm—ls 


2 
tn COSVY OK 2 PÜAm-—ı + 0 %m + B) cos v j om—ı 


RU Er Ne 
Mer ee om —ı + 0m + PB) sinv jom—ıl. 


Mit den Bezeichnungen (2) ist 


a, (ß) 1 a, M) p—1 a, (0a, +B) h b, ("a +P) r 
e h = D ei h 7 2 ei B I a „ om] 68), 


ähnlich 
b, (B) 1 b, \Dd o—1 b, (ca, +B) a (am +B) 


Mit einer positiven ganzen Zahl m sei 


TB) (,,) os |.) sin’». f 72 er en N E | 
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Dann erhält man aus (5) 


7) \(e%) 
Ina) =( ; a ve| x - sinvoo, 


p" s 


ri a, (09) »—1 a, (0%, +5 39) Mr 
 p \\pri1s co8y day, + 2 Ki ) COS 9 0 Om 0089 00% 


b, Memtee) 
> Dez Ss 


2»—1 b, (va, +9 %g) ! ii a, \(oo,+0 %g) ! TH £ 
— 5 mis c0Sv 0%, SINP OA, + sın v0 Sn v0 &, 
1 


j e Due > 
sın v 0 Om, cosvoa,| — pm—1g smv 0a, 
Paz s 


er! a, \e sg) a BoAd Sa) Se 
ln, 089 0%, — pm—ı ) sinvoo, 


Pl (01 +9a,) (Oam+oa 
+2 je ln ke sta 


Mit der Bezeichnung (6) ist 


1,0%) = 


sn Pa S 
ee &) + > Tm—ı(0 m + © ee) ee a (2): 
Mit Hilfe von (7) findet man 


| 5 )= ae Me 3 Tre) 


p"s pP DE oe 
Ii—1 p—1l»p—1 »—1 p—1 
WRITER 2 2 > Seue), Tal &n + 
P r=1 (dir da +.) a=l 1 ,-1 “Reg, 
EB 
ne er ee ie 


wo IQ <I!I<n) eine ganze und positive Zahl bedeutet. Die Bezeichnung (i,, iy, ..., L,) 
e{2,3,..., 2} bedeutet, daß die Indizes {i,,i,,...,i,} die ( ) Kombinationen von r In- 


dizes, die man mit den Zahlen 2,3,..., ! bilden kann, durchläuft. 


Um die Formel (8) zu beweisen nimmt man an, daß sie für 1 = I befriedigt ist. Man ersetzt I 
durch ! in (8) und man formt die Ausdrücke -T,„_, mit Hilfe von (7) um. Man erhält: 


| a, ) Bi A ,) de SE Ta—1-1l0,) + —. > 5 Te ee) 


pP" s P pt! a=ı ptia=ı = 
1 il »—1 p—1 »—1 py—1 
+ a oa £37 > > 3 Br PR 1. 72407 677 + 
Des 1 (ds ber, ß,) a=1l 1 5,1 o=1 
ei2,3,...,7} 
+ 9% nt + mitt + 0 in) 
1 !-1 Pet 91 »—i p—i 
Sr 2. 5 Ta—1-1(%7,1% PN] Tr 
p EUeZ1 (dr dar 2 23,) a1 5,1 1 9, „= A 
e{2,8, ...,1} 
ee ar N 0) en Sara sl). 
Man wird weiter die symbolische Gleichung: 
edel} 
Benin... tree 0,30..0 #0) 
benützen. Mit dem letzten Glied werden die Kombinationen von r Indices bezeichnet, die man 
bekommt, wenn einer von ihnen den Wert ! + 1 annimmt und die Kombination fi, ig... .,1,—ı} 
die | Kombination von r— 1 Indizes durchläuft, die man mit den Zahlen 2,3,... Ai 


Bl 
32* 


2 * 2 
er (97 ce FO: el 


u ee a=1 “m De z ai 7 Bis br a 
ü «2,3. 5, 2 = i 


2° IR 


i 
ae Fr #2 (rien N er >, Er 


e{2,8,. ul 
+ + nt te, IR 3 
Wenn man r=_+ 1 setzt, wird die vorige Se: zu: 2 
(el len) to 2 ee m Turlamntenan 


£ u 7 
1 1 2—1 »—19-19- 


1 L g j 
= 2 E” E 5 5 :-Taoi-ıt 0, tn 
pri (via...) 1-1 09-1 0,-1 05-1 Zn, 


EOREN SE SE 
za p—1p—1lp—1pr—1 9—1 - u 


ne Tail + 


+1 oe=1 (ü, ee ip +1) a=1 a! 1 Sa jr 
ei2,3,.. ah} 3 « ur ci ® 
rm tr mit Hrn) - 


Wir haben die Gleichung (9) bekommen. Wenn die een: (8) für Z = I richtig ist, ist 
sie also auch richtig für = 1 +1. 
Für !=2 wird (8) zu 


(| A, +3 az, Tan) + 5. 2 Tan +, m-ı)- 


p"s Balz a Sr 


Aus (5) ergibt sich für # = 0 und m = n mit den Bezeichnungen (6) und mit Hilfe von (7) 


ar 1 da, 1 a, 129-2 
® a = Ende +7 7 Di T2—ı(0ı %y) en en fi tz u Ta—2(51 &) 
1 vr—ıipr—i 
er 2 Tn—2(01 1 + In). 


p» o=1l gel 


Damit ist die Richtigkeit der Gleichung (8) bewiesen. 
Wenn man in (8) = n einführt, bekommt man: 

a, 1 a, 1 es e n—ı p—1lp—ıip—ı p—1 
re) tz Tantra: 2m 
P’s pP S pP u Pe (Ur iarı..,,) a=l 5 =1o=1 a \ 

ers a z ' j 
+ 0% +10, m—iu+ı tt: + Oi, n—i.+1) a Ba 2 ee SE 
Aus (6): 


T,(ß) = WE cosy ß — RR kin va: valsikbeni er RN | | | 
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Wir bezeichnen: 


t == =E p—1lp—1p»-—1 vr—1 
On = m ei ee 2 Zu oe Z Toon + 
in 


r=1 (dir ia» 3 *,8,) o=!1 5-1 ,=1 


E{2,3,...,n} 
+ 9, n—iu+1 + 0, n—i,+1 +. + Qi, | ee Me en (13). 
Damit geht (11) über in 
een 14 
BE an 
und wir erhalten endgültig 
a, 1 ja, en 
ea Z ET A (15). 


Für n = 1 hat der Ausdruck (13) für 0, keinen Sinn mehr. Um 6, zu finden, führen wir in (5) 
n=1lundß=0 ein. Mit der Bezeichnung (12) erhält man 


WEHEN 


Ps P a=1 
Daraus und mit Hilfe von (14) bekommt man 
1 
0, were: ei T (0, 9) 


- a,\® a, /b,\®) 
Die Gleichung (15) ergibt den gesuchten Ausdruck von als Funktion von 5) ; Ei > 
wobei ß die durch (13) und (15) gegebenen Werte annimmt. P 


In ähnlicher Weise beweist man 


Br] 
Br 


1 
On = Zn z So(01 &) 


n—1l p—-ip—1ip—1ip—1l 
PH = De zZ So (Hm + 0, nn +1 + 0, m Hı tt ee] (17). 
r=1 (ir in - Pr „) 11 9, at 9, =1 0; =1 - 
e12,3,...,n} ü 


(Für n = 1 verschwindet das Glied von o,, das die vielfachen Summen enthält.) Außerdem ist 


sw =(F)" eos»ß +(*) sinvA KEITEN: 


Wir werden jetzt noch andere Ausdrücke für ei und (7) aus (15) und (16) herleiten. 


Wir betrachten die Größen a; a,» 25... 5@äsAys* ++ 529» WO die Indexkombinationen 


BDA Ale en al Ayo Ass pre An) die 1) (2): EN, (n) Kombinationen von 1,2,...,n Indices 
durchlaufen, die man mit den Zahlen 1, 2,...,n bilden kann. Die Größen qa,,,.,...,a, haben 


definitionsgemäß die Eigenschaft, daß sie sich nicht ändern, wenn bei ihnen die Indices A,, As, 
.,A, vertauscht werden. Man hat 


n dl N 
2 (« ei 2 Ag, 5 -u+l,0—istl, ...4 4) =„2 & EN e 19). 


o=1 r=1 (dr dar «+ +3.) o=1 (Aı A 2) 
ED Byrne zL0h ENTF. n} 
Für o =1 wird das Glied 
o—1 
Lo 3 Mo-üut+lo—utl...,o—i+1 
r=1 (ds ie en Ä. 4 


e{2,3;...,0 
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aus (19) weggelassen. Um Gleichung (19) zu beweisen, bezeichnet man das erste Glied mit D: 


n—1l n 
D= 504 Pa Do A,o—ht+l,o—itl,...,a—iptl: 
o=1 r=1l o=r+1 ER i 
et, 3,. ai 
Die Summe 
DO ar o— ht... Hl (20) 
o=r+1 (ivine..ni,) 
Eei2,8,2..,0) 
hat 
I 
BU. r ® ar rtr1 
Glieder. 


Wenn o, < o, ist, unterscheidet sich die Indexkombination (iur. A-ir+1) von 
der Kombination (, a» -ı +, ga —haht+t1l..,%—J+ 1); hier sind (I, in..., ir) und - 
(for +. /r) Zwei mit- den Indices 29,4, io und 2, 3,9. gebildete Kombinationen. 
In der Tat sind die Indices 0,0 —ı +1,...,0, —i,+ 1 alle kleiner als Oz. Wenn die mit 
den Indices 2,3,.... gebildeten Kombinationen (i,iy,...,i,) und (ji ]» - .. j,) verschieden 
sind, sind (ao — u +1...,0o—L+1) und („o—ı +1,...‚,0—j,+1) auch ver- 


schieden. Also ergeben, da die Glieder der Summe (20) aus ( untereinander verschiedenen 


” 
r+1 
Kombinationen von r + 1 Indices bestehen, die den Zahlen 1,2,...,n entnommen sind, alle 
Kombinationen, die aus r + 1 den Zahlen 1,2,...,n entnommenen Indices gebildet werden, 
die Summe 


au} 799 PERERDIPe VOR 
(Auda set) 
EL N) 
Mit@=r-1 erhält man 
hi wi 
D= >>; OR a Be 
Heihdh...,m P=2 (huda...,A,) 

e{l,2,...,n} 
wo E das zweite Glied der Gleichung (19) bezeichnet. Damit ist die a, (19) bewiesen. 
Wenn man den Indices A, Ay ...,4, "Werte aus der Reihe der Zahlen 1,2,...,n gibt, ist der 
Ausdruck: 


p—lp—1lp—1 
ON a a Dre Tel 0 + 420, ++ 4,02,) ee LE 


m=l m=1 u,=1 


eindeutig bestimmt, und a, ,.,...., 32 verändert sich nicht, wenn die Indices ra ER ‘ unter- 
einander vertauscht werden. Tucß) ist durch (12) bestimmt. 
Wir setzen jetzt 


(deg) D),\0°0) 
2 A, : R De 
At 0-(®) os”) sin vr ai AITTERAT TREE 


“ b,\0°0) h de). 
BV*°ı) -(*) cosy ja, + (*) SINiV.] gt un hr nn EN 


Dann ist aus (12) und (12a) zu ersehen, daß AP’ — T,( &g) ist, und es folgt: 


a 04, haar Hp °2,,) 2 


= Tom ++ Mo 0%) ee 


Aus (13) und (15) erhält man, wenn man o durch 9 ersetzt 


r=1 (Ur ia...) o=1l 1 %,—1 
RC 


A, ıl N Sm n 0=1 pn 
Ei le +, a „2 T (0, &) Tr 2 > 3 PIRRN en En 65% (T,1 0 + 


TR + + On ar + Oi, | . 
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Mit den Bezeichnungen 21) und mit m = 0, 4 = Oi, .., Urt = 0, ist 


da, il, o—1 
Kers, = | +2 2 ‚(@ 2 PA j en) e 


p" r=1 (dis dar - 3%.) 


e{2,3,...,0 
Mit Hilfe von (19) und (22) und mit o = 9 ergibt sich 


a, lg: = De 
En a a ea 
pP” Ss pP S Pl (Ayvdmı..d,) al mel Hg 

€e{2,3,...,n} 


Ersetzt man 4» .- -» 4, durch Op Mr On so bekommt man 


Sala Eee 


n n 
[u p o=1 (Av Aa... Ay) o,-1 u 


Auf gleiche Weise kann man die Ausdrücke für en und Ei finden. Man findet schließlich: 


EN Ay er a N 
I ,) eh nn (hg |. 


a 1 A,\ a Dil o9.0+072 04, +:.:':-+0, & 
el aa Zu A Ay Ay) 
3 e h p” | x + 2 a BR nr u he ee (23). 
b 1 b,\ N, p—1 0,04, +0,,0,,+ "+0, a 
re BOCH Art Lu) 
Pe ,) p" ( he 2 en er Sr 


In diesen Gleichungen sind A und B durch die Ausdrücke (12a) und (12b) gegeben. 
Für n = 3 bedeuten die BE BE 


FRE (+2 > AR I Aa 4 5 Aa + a Augen 


3 ms 
p 3) pP ah=1l 0-1 0-1 09,03 =1 


ze 
+ Ben A(92%2+03 83) + u A'd%ss+0n aı) + zz Acı aı +03 03403 -j - 


0, 0=1 0,0=1 01, 03, 051 


a=1 0-1 o=1 0,01 


i b, 1 b, Da v»—1 p»—1 »—1 
| 3 = | )+ > Bias) + Pa Base.) 4 er Base) 4 > B(6ı 1 +02 02) 


p—1 y—1 »—1 
+ 7 B92 2403 0) En > B(03 as+01 61) + RR B(cı 1 +02 Ba! 4 


03,051 0,0=1 01, 09, 09 =1 
Der Gleichung (3) entnimmt man: 


27 27 
> 2% = ng at DRTÄ gR (An) er 2 
Mit Hilfe dieser Gleichung kann man die Werte von 4, %, .. . ‚07, aus (23) als Funktion 
der kleinsten von ihnen, derjenigen also, die den größten Index besitzt, ausdrücken. 


Auf Grund der Formeln (23) kann man die Näherungswerte der Koeffizienten t 


a b er ee i 
| 2“ - % bestimmen, die einer Einteilung des Periodenintervalls in p*s gleiche Teile ent- 
DESIINDES 
i . N a a b . 2 EEE 
sprechen, wenn man die Werte der Koeffizienten ( ‘| R | Äh | ) kennt, die zu einer Einteilung 
s s 

des Periodenintervalls in s gleiche Teile gehören und man von verschiedenen wohlbestimmten 
Anfangspunkten O0 auf der Abszissenachse O u ausgeht. 


Ay 


b - i 
Die Bestimmung der Koeffizienten () (2): 3 für s = 24 gleiche Teile kann mit 


‘den Tabellen von ZıprzR£r durchgeführt werden. 


Je größer die Zahl der Teile p*s ist, desto genauer sind die Werte, die man für die Koeffi- 


zienten ( ° \, ( ” ) | one (23) erhält. 


Die Vorteile der Methode der Bestimmung der angenäherten Werte der FouRIER-Koeffi- 
zienten, die in vorliegender Arbeit beschrieben ist, sind die folgenden: 

Die Berechnung der Werte der Koeffizienten der Fourer-Reihe mit den vorhandenen 
Rechnungsschemen (Runge, ZIPPERER) benötigt sehr ausgedehnte Rechnungen, wenn eine 
größere Genauigkeit für die Werte der Koeffizienten verlangt wird und zu diesem Zwecke die 
Periode in eine größere Zahl gleicher Teile geteilt wird. 

Mit Hilfe der Gleichungen (23) kann man jetzt aber die Werte der FourtIeEr-Koeffizienten 
mit einer großen Genauigkeit erhalten, wenn man von den nur ungenauen Näherungswerten 
dieser Koeffizienten ausgeht, die man sehr schnell und einfach erhält, wenn man eins der vor- 
vorhandenen Rechnungsschemen (Runge, ZIPPERER), für eine Einteilung der Periode in eine 
kleine Zahl gleicher Teile anwendet. > 

Die Gleichungen (23) ermöglichen die Berechnung der angenäherten Werte der Koeffi- 
zienten der Fourter-Reihe für eine viel größere Skala der Teilzahlen der Periode als die gewöhn- 
lichen Rechnungsschemen (RuNGE, ZIPPERER), die für einen bestimmten Wert der Teilzahlen 
ausgearbeitet sind. B 

Diese Gleichungen ermöglichen sogar die Berechnung von Näherungswerten für die FOURIER- 
Koeffizienten mit einem beliebigen Annäherungsgrad. 

Der Weg den man einzuschlagen hat, um die Fovrıer-Koeffizienten a,, a,, b,, einer periodi- 
schen Funktion mit Hilfe der Gleichungen (23) zu berechnen, ist demnach der folgende: 

Man wählt die Zahl s der gleichen Teile, in die man die Periode zerlegt so, daß man schnell 
und einfach die genäherten Werte der Koeffizienten a,, a,, b,, mit den Schablonen von ZIPPERER 
oder mit den Tabellen von Runge rechnen kann. Danach wählt man die Werte von p und n 
gemäß der gewünschten feineren Unterteilung der Periode in p"s gleiche Teile und rechnet die 


Werte der Koeffizienten ® )- 5 ) | o ) mit Hilfe der Gleichungen (23) aus. Dazu werden 
s Des 1 

die Werte von a, mit Hilfe der Gleichungen (24) und für jedes x, und für den angenommenen 
Wert von s die Näherungswerte der Koeffizienten a,, a,, b, für den Fall bestimmt, daß man 
der Periode O0 bis 2% eine durch die Gleichung (24) für &, bestimmte Translation entlang der 
Abszissenachse O u erteilt. Die gefundenen Werte von &, und der Koeffizienten a,, a,, b,, setzt 
man in die Gleichungen (23) ein und erhält die genaueren Werte dieser Koeffizienten a,, a,, b, für 
eine Zerlegung der Periode in p” s gleiche Teile. 

Man kann auf diese Weise die Koeffizienten der FOURIER-Reihe einer periodischen Funktion 
für jede Ordnungszahl » berechnen. 


Bei der Rechnung kann man die Endgleichungen sowohl in der Form (23), wie auch in 
der Form (15) und (16) verwenden. 
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Large Strain Creep of Rotating Cylinders”) 
By F.P. J. Rımkorr and J. R. Luxe 


In der vorliegenden Abhandlung werden Gleichungen für die Kriechgeschwindigkeiten, die Kriechdeh- 
nungen und die Spannungen in rotierenden runden Hohlzylindern aus isotropem und homogenem Werkstoff 
unter folgenden Annahmen hergeleitet: 1. Die Wichte des Werkstoffes bleibt unverändert. 2. Die Verhältnisse 
der Hauptschubspannungen zu den entsprechenden Hauptschubdehnungsgeschwindigkeiten sind untereinander 
gleich. 3. Die Vergleichsdehnungsgeschwindigkeit ist eine Funktion der Vergleichsspannung. Um das Ver- 
halten :des Zylinders auch bei größeren Verformungen zu erfassen, wird ein logarithmisches Dehnungsgesetz 
verwendet. Auf diese Weise kann die Zeit bis zum Versagen des Zylinders berechnet werden. Den Gleichungen 


liegt das Kriechgesetz & = B 0" zugrunde. 


Using the generally accepted assumptions of constant density, of equality of the ratios of principal shear 
stresses to the corresponding principal shear strain rates, and of existance of a relation between the significant 
-strain rate and the significant stress, equations are developed for the strain rates, strains, and siresses of a 
rotating hollow eylinder made of an isotropic and homogeneous material which is subject to creep. T'he strains 
are considered large, requiring the use of finite-strain theory. T'he application of the finite-strain theory also 
permits the prediction of the time at which the cylinder will fail. The creep law € = B 0" is used in the 


derivation of the.equations. 


B cTaTbe BbIBONATCA YPaBHeHHA HI CKOPOCTU YTEeYKH, PACTAReHUA B CHJIY NOA3ZYUeCTu 
Hu HANPSBREHNÜ BO BPAIMalI0UXCA IIOJIBIX KPYTOBbIX NHJINHHAPAX 43 U30TPOIHOTO U ONHOPONHOTO 
marepnasa. Ilpı 3ToMm npenmonaraetca: 1) YıeibHbIä Bec MaTepuasıa He MeHsieTch. 2) OTHO- 
IIeHNHA HAIPsBKeHHÄ TIIAOHOTO CABUTA KR COOTBETCTBYIINHM CKOPOCTHM YTEYKU PABHbI M&KNY 


co6oü. 3) CROPOCTB HAUPSBREHHA CPABHEHUA eCTb PYHRIHA HAIPAHKeHUA CPaBHeHHuf. 


UTO6bI ONNCATBL HOBeNEHUA NUHAINUHNPA Taryke pm 66NbBIIMX MeDOPManuuAX, UCHOABZYETcH 
JoTapuhMHyeckmü 3aKOH pactazkeHun. TakuM 06Pa30oM MOJKET ÖBITL BbIYHCJIEHO BPEeMA IIPe- 


Kpammeunna neictsun umamanpa. B 0ocHoBe ypaBHeHnüä Je&KUT 3akoH TeryuecTm & = B 0". 


Nomenelature 
Symbol; 
In Natural logarithm (to base “e’”) 
E True (natural) strain 
€ True (natural) strain rate, also significant true strain rate (hr=2) 
o True stress, also significant true stress (1b in?) 
K Modulus of creep (lb in? hr) 
B Creep coefficient (in? 1b? hr-1) 
n Creep exponent 
r Radius before deformation (in) 
7% Radius after deformation (in) 
7 Change in radius due to deformation (in) 
2 Time rate of radius change (in hr!) 
a Original inner radius (in) 
b Original outer radius (in) 
a’ Inner radius during deformation (in) 
b’ Outer radius during deformation (in) 
F Force (1b) 
H Ratio of the significant true strain rate to the axial true strain rate 
& x evaluated at the inner radius 
ß x evaluated at the outer radius 
R Ratio of the inner and outer radii during deformation 
Ro Ratio of the initial inner and outer radii 
t Time (BrRe 
0 Mass density (lb sec? in-*) 
10) Angular velocity (sec): 
L Bow? (b? — a?) (hr!) 
Subseripts 
[0 tangential 
5 radial 
2 axial 
[2 at inner radius 
b at outer radius 
f at failure 


1. Introduetion 


In problems concerned with deformations in the ranges of elasticity and primary creep 
it is generally sufficiently accurate to make the assumptions that strains are small and that 
they can be referred to the original dimensions of the body rather than to the instantaneous, 


*) Auszug aus der Magisterarbeit von J. R. Lukr, Toronto 1960. Ein Auszug aus dieser Abhandlung 
wurde auf dem X. Internationalen Kongreß für Angewandte Mechanik vorgetragen unter dem Titel ‚On the 


Creep Behaviour of Rotating Cylinders“. 
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deformed, dimensions. During secondary and tertiary creep, however, the magnitudes of de- 
formations may reach values where the assumption that strains are small can no longer be made 
and the application of the finite strain concept becomes mandatory. 

Based upon the large strain theory, an expression is obtained which relates the creep rate 
of the cylinder to the creep deformation and the angular veloeity of the cylinder. Integration 
of this equation leads to a relation giving the creep deformation of the cylinder at any time for 

a given angular velocity. t BR 

| The creep failure time is then obtained by letting the creep deformation approach infinity. 
It so happens that the time necessary to reach infinite deformation, turns out to be a finite 
quantity. The time to reach infinite deformation, called creep failure time, can be used to prediet 
the time to actual fraeture. Naturally, the cylinder will fracture before it has reached infinite 
deformation. The expansion of the cylinder, however, proceeds so slowly at the beginning and so 
rapidly towards the end of the life of the vessel, that the time to fracture and the time to reach 
infinite deformation differ only by an insignificant amount. This means that the creep failure 
time can be used to predict the time at which the cylinder will fracture. 

The creep failure time for the case of simple tensile creep was predicted theoretically by 
Orowan [1] and Horr [2] and shown by the latter to be in good agreement with experimental 
findings. Creep failure times for the hollow closed-end eylinder and the hollow sphere under 
internal pressure were obtained by Rımrorr [3] and Rımrorr, Mırrs and Marıs [4]. The present 
paper extends the theory to the problem of the rotating cylinder. 


2. Assumptions 


In the stress analysis of elastic systems in equilibrium, the linear relationship between 
stress and strain permits the use of the method of superposition. Under creep conditions however, 
the nonlinear relationship between stress and strain precludes superposition and the stresses 
must be considered as a combined effect. The following is a compilation of the assumptions and 
definitions on which the present study is based. 

1) The cylinder is made of an isotropie and homogeneous material and retains its circular 
shape during deformation. 


2) The volume of the material remains constant. Volume constaney can be expressed as 


sts. +3 =0 .. 0, Br Br (1). 
3) The true axial strain is constant with respect to the radius. That is 
ge (2). 


Imposing this condition limits the discussion to long eylinders. 


4) The principal strain axes do not rotate during creep. This permits differentiation of 
Equations (1) and (2) with respect to time to give the following ereep strain rate relations 


tee (3), 
ee (4). 


5) The ratios of the principal shear strain rates to the principal shear stresses are constant 
at a given point in the cylinder but may vary from point to point and with time, 


Ha u 3 (5) 
in rr “ = = x “ > . * . . * ” 
9% —-% %—-% 0, 2K 


Equations (5) are a statement to the effect that Monr’s cireles of stress and strain rate are geome- 
trically similar. 


6) There exists a significant true strain rate versus significant true stress relationship which 
can be obtained from a simple tension test. That is 


&=flo) 2»... . 5 000 MOSE BE 
In this paper, the above relationship will be taken as 
ET 


7) The significant true stress will be given by 


1 
= 7 [0-0 + (a a ee 


Equation (8), the von Mises stress invariant, is felt to give the best correlation between the 
multiaxial stress system and the simple tension stress system for creep. 


£ 
E 
R 
2 
- 
f 
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8) The significant true strain rate corresponding to the above significant true stress is 
given by 


De a N (9). 


3. Compatibility Condition 
In the following work, the development of a number of equations based upon compatibility 
of strains is given. These equations will be needed later in the analysis. 
The true circumferential strain is defined 
as (Fig. 1) 


&,—= In ( +7)- In(% 
r 


Fi 


) E10: 
The true radial strain is defined as 
du dr’ 
nl). ih 


Equations (10) and (11) may be differentiated 
with respect to time to give the true strain rates. 


The true eircumferential strain rate is then 


— er ur RE FRE BF 1 
‚-T (13) 
Combining Equations (12) and (13) and BEFORE DEFORMATION | AFTER DEFORMATIOR 
eliminating r’ and dr’ gives the compatibility 
equation 
r’ de A ; 
dr’ ’— Er &, URN „ie (14). Fig. 1. Cross Section of a Cylinder 


Substituting Equation (4) into Equation (3) and solving for &, gives 
BB el De (15). 
Substituting Equation (15) for &, into Equation (14) gives 
L. de, 
dr’ 


Equation (16) can be integrated by separating the variables and results in 


R lade : 
- A EEE ee 17), 


where c is a constant of integration with respect to radius but a variable with respect to time. 
Solving Equation (17) for c gives 


ee A ed > AD Ta re (18). 
For r’ = a’ and r’ = b’ Equation (18) becomes 
BO EEE BEER) er. een (19). 
Substituting &, from Equation (17) into Equation (15) gives 
& INFe . 
u, (ath) Be co Et: (20) 
Substituting Equation (15) into Equation (9) gives 
3 Br an ; 
>= \5 En un. u un. no (21). 


“ PN Aa BEINE 


A ie } 

r TE Te 

Say DB IRA Ve 
\ 


4.08 
u e analysis a 
in eonjunction with rations pı | 

N equilibrium equation for aro ating cylinde can IOW 


= € 


r do, lan r ThR 
ı “ e rn a - en Mu * Hr kur Bed ee 
Ur et .. 
. gateriah ER 4 r a 2 r 5 
Equation em) en rearranged Dora « (£r) Ima (EI) rei a 


pr hen Ruh“ ige DARF De 
(0) gar a 


DT 


R r = r vi vu D; & 
Combining Equations (3), (4) and (5) gives En r Kr. 
1 In; on SIR Dr 
= +0) Er: * j 
1 1 er 
= TEE I 
ze 2 oe 
u= [950% + 2 ee ae ee . a) 
where rain - 
o | 
K = = Bm 0. 
Solving Equation (31) for o, and substituting from me: (32) gives the Inga ‘ 
stress j KW 
= k+,@,+0) tunen de en 2.8). 


Substituting o, from Equation (33) into AN (8) and solving for (0, — 0,), 


(0%, —0,) Smras kayur Tr ee Be; (34). 
By rearranging Equation (16) it is possible to obtain 
dr’ , —de, 
(2:,+k) 
Solving Equation (21) for &, gives . | Bu 


ee... a (36). 


i 
E 
3 
H 
Br, 
# 
a 
= 
= 
a 
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Differentiating Equation (36) and recalling that kis a function of time alone, 
ae Eee 
2 (22 — ka)u2 
Substituting de, from Equation (37) and (23, + k) from Equation (22) into Equation 
(35) leads to 
dr’ : & de 38 
= Dita see, 


Returning to Equation (28) and substituting (o,— 0,) from Equation (34) and = from 
Equation (38), 


— ode MT 
do, = Ber Te a 0 (0% 'e dr or ee We ee (39). 
y3 (& ei; kayu2 


Equation (7) may be rewritten to read 


Sp 
0=|5 RENT a ae Bhaaı )> 
Insertion of the above value for o into Equation (39) gives 
Be START: 
do, = — e 2 23 BIDNPHLTEENUDH 20 DIENTE NN, 
V3 Biln (& — RyUl2 z 


Imposing the boundary condition that o, = O0 at the inner radius of a hollow cylinder and 


then integrating Equation (41) gives 


& 
—1 ellm de 2 
iR ie An Beten 3 BEN ee ee a 
Y3 Bin (& Au kayu2 2 
Ea 
When the upper integration limit is &, or b’ then o, is again equal to zero and 
& FE . 
elln de 3 
== 1/n Hide 
Ir: 5 Bin 0 @2 (b eh iz er A: 


€a 


5. The Axial Force 
Another independent equation involving functions of & and r’ is required for the solution 
of the problem and this is obtained by imposing the boundary condition that the resultant axial 
force acting on the cylinder be zero. 
Returning to Equation (28), substituting (o, — 0o,) from Equatiom (34) and integrating, 
r’ 
2.2002 a We, r 
a | (ey > Le ee ee ıry 
a’ 


Rearranging Equation (34) gives 


+0) - ne K2)U2 + 0, IE: B 2. uldn); 


Substituting o, from Equation (44) into Equation (45) and then substituting the resulting 


| equation for (o, + 0,)/2 into Equation (33) gives 


r’ 


N . dr= — 0.09? 
=, Basıfref 2\1/2 ar E77 Be. Ja\ll2 en RO NR Be: 16). 
G..= kt (@— — k) +7 |2e K?) r' D) (r @2) (46) 
Since the resultant axial force is zero, 
b’ 
Bde 2. u. sn. a. (di 
a’ 


Thus n 
ne N... Ne. AB). 
a’ 
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Substituting Equation (46) for o, into Equation (48) gives 
vr b’ 9 b r dr’ 
ee a I (@ — ey —| rar 
- [Feier te 
a a a’ a f 
2 2 
23 = (rad) Te. ee er a 4: (49). 
a’ 
The double integral can be integrated by parts, giving 
b’ 
klan n se ig _ 3 SE (Brad 2.2.2.2... (50). 
Ei 
a’ 


Integrating Equation (44) between limits a’ and b’ makes o, equal to zero, hence 
b’ 
Im a ER 
—I1—- (a2 — RR = (fa .. 2. 02 1). 
„je in 2 aa 61) 
Substitution of Equation (51) into the SIE for the double integral given in (50) gives 


2 
2 b'2 (r2 — a2) — — JE Se ZA (52). 


By evaluating the last integral in Equation (49), substituting (52) for the double integral 


and simplifying, one arrives at 
b’ 


# 2 
= krdat=— ru (1 ee (53). 


From rearranging Equation (25), 


z c 
I 3(@_— gr ee - (54). 
From Equation (38) 
drug 1 & de 
= ae) a a a nn DE 
Thus 
it a 


2 y3 (&® — ke) 
Substituting Equations (55) and (40) into Equation (53) and rearranging leads to 


Ep 
an de 1/n 
VE: LOL 
Ea 
Considering Equations (23) and (24) it can be shown that 
b’2 + a’? 1 1 1 
| e )= 7 (m Eee ER] N 


‚2 


- Bar u 
Substituting ( -) from Equation (57) and & Bin 9 w® (b'? — a’2) from Equation (43) 
into Equation (56) gives 


&) PR 


auf, ein de 1 1 zum dä 
— 2y3k i =” 
/ kn — k2)312 ut 2 — jayu2 "a ee ED 


Eu Eu 
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R At this point it becomes advantageous to define and introduce the following dimensionless 
ratıos 


Eu 
Ne ee a ae ee en u, 
h (60), 


ß -2 EEE ARTE (61). 


Using these new expressions, Equation (58) becomes 


B B 
I im da 1 ıl yılm dx 
—- —— eng) ee er — u EEE 62). 
u me) (21 = 


From rearranging Equation (10) 


a a ee ER (63) 
Thus 
ee Pan 9 ES ee (64), 
De ae ER (65). 
From considering Equations (1), (2), (10) and (11) it can be shown that 
ET N ie) er en: (66). 
Thus 
a (ie a A er (67). 


Dividing both sides of Equation (67) by b’? and employing Equation (65) gives 


28 
ER En 2 en Res (68). 
I) 
At this point, let us define and introduce the dimensionsless ratios 
a’ 
EEE VEREINE, 69), 
r (69) 
Bi rn EEE BEE N (70) 
Equation (68) then becomes 
1—R 
Pearl re ae ed 71 
Similarly, it can be shown that 
‚, 1—-R N ö 
28, ,+% 
EL R IR Sue a (72). 
Differentiating Equation (71) with respect to time gives 
. 2, Fed SER 
9 — . RE WI 


Equation (24) can be extended to include the new quantities R, a and $ and becomes 


RE el 7} 


are 


DD“ 92 Eu tk (&2 — keyu2 = Vo® = 


| WaThe ereep x 

i ea Aveo re 

E reasons which will become evident later, it ) Sc 
r with a particular solution. In the case at hand, n 


2 is A than 1 will be inves ed later. 
FUN T = 
6.Sotinernel 0000 


When n = 1, the integrals in Equation (62) become PR 


dr A % 
Ye—1 - Ge rot rübirikl u 


& 


Making these ne into Equation (62), replacing /ß? ve R?ya®— land soling 


ee ‘ 
for Ya? — 1 gives £ 


a = 

a a ee a a NE (80) 
and ö 

—  —2y3R? 1 Acta a j 

a a a de ee a (sl). 


The relationships expressed by Equations (80) and (81) are plotted in Figure 2. 
Considering Equations (75) and (80), Equation (77) becomes ö 
—1(1+R) dR . 


ln a 
3AL— BERN, SR a re A; 
Integrating Equation (82) gives 2) Bein 
ne iR 
PR en Va 3 
ss Vame nn (83), 


from which e 
% R 1 . R2 148 
k— Me 84). 
. je ZRNTR = 
Considering Equations En and (67), Equation (43) may be written as 


kım * m dy 


Teer = u! Bin oa (a Bd) ERTL TR EEE (85). 


x 


In this case, n = I and Equation (85) can be integrated using Equation (79) to give 


ERURZI- Vai) = Bora NUN (86). ) 
| 
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Substituting e* from Equation (84), k from Equation (82), YP— 1 from Equation (81), 
Ya® — 1 from Equation (80) and solving for the time {by on gives 


Ah 4 R, \8 [ /1 — R2\u8 
ee j ar 2. en). 
i 2:7 
Evaluating Equation (87) by series integration and making the substitution 
gives Li DB 0 Dee (88) 
1 Be 64—1 
N ge 
Be |" 5 AH: ) 
‚ Fi 1 ei. in je 1 El (89). 
(A)! yo I(63— 1) 
R, = a/b and R=a')' 
o 0.1 0.5 1.0 


L« Bpu2(lh? -a2) 


R, = a/b 


Fig. 2. Values of Ve — 1, Yß®—1, o, and ß versus Rforrn=1 Fig. 3. Creep Failure Time versus Initial Radii Ratio forrn=1 


Equation (89) states that if the initial radii ratio R, and the constant L are assigned values, 
the instantaneous radii ratio R can be determined at any time in the future. With R known, 
the cylinder can be completely analysed by employing equations developed previously. 

Mathematically the creep failure time £,, is the time at which the creep strains reach infinity. 
Physically, failure by fracture will occur when the strains are much, much less than infinity. 
However, because the strains change extremely rapidly from relatively small values to infinity 
at times close to the theoretical failure time, this theoretical failure time f; becomes a useful 
approximation for the life of the cylinder. In the case at hand, the condition that the strains 
be infinity, is equivalent to the condition that the instantaneous radii ratio R be equal to one. 
This more convenient latter condition will be employed. 

Thus, from Equation (89) it can be seen that 


„sothter) 


4=0 a3) Her 1) 


Values of the dimensionless creep failure time parameter L - t, are plotted in Figure 3 for 
cylinders having an initial radii ratio R, between 0.1 and 1.0. With the value of R,= a/b known, 
33 


Legale — pi . (90). 


TIME PARAMETER L*t 


Der 
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Figure 3 will provide Lt; and sinceL is a constant known from initial conditions the creep failure 
time 1; can be determined and the life expectancy of the cylinder found. 

In the solution of the problem, it became apparent that finding the stresses and strains as 
functions of time directly would be much more complicated than employing a two step procedure, 
the first step of which involves the determination of Rat a given time and the second involves 
the determination of the stresses and strains as functions of the now known quantity Fr 

For a given R,, Equation (89) allows the determination of R as a function of the time 
parameter Li. As an example, when R, = 1/4 
Equation (89) becomes 


ol) 


1.4 
Deu. 


L = Bpu°lp2 - a2) 


L-t = 2.58 1 
abs a3 —)162— 

\ (91). 
Values of L - £ obtained for the range 


0.25 < R< 1.0 are plotted in Figure 4 and ta- 
bulated in Table I. 


1.0 


TableI 
1 
1 T 
0. R L t Epa Era k 
0.25 0 0 0 0 
30 | 0.215 0.2312 —0.1608 —0.0704 
.40 .566 .6215 — „4285 — .1930 
‚50 .810 .9580 — ‚6530 — .3050 
.60 .980 1.275 — .8550 — ‚4200 
.70 1.120 1.600 — 1.055 — ‚5450 
.80 1.240 1.990 — 1.284 — ‚7060 
.90 1.280 2.560 — 1.600 — .9600 ö 
.95 1.305 3.045 —1.855 —1.19 
A 0. .96 1.310 3.210 — 1.940 —1.27 
on .97 1.315 3.415 —2.045 —1.37 
INSTANTANEOUS RADII RATIO R = a'/d* .98 1.320 3.700 —2.190 —1.51 
.99 1.325 4.170 —2.430 | —1.74 
Fig. 4. Instantaneous Radii Ratio versus Time for R, = 0,25 1.00 1.330 e) —'&9 — © 


Strains and deformations can be found as functions of the ratio R (and hence as function 
of the time 2). 


The axial strain is found from Equation (83) as 


ER R, 1— R? 
Iren rt "R en 
From Equation (84) 
RB, . ir Raun 
Krss Tue 
F ee: BEN 
From Equation (72) 
„2 artk c = R) 29 
e% al a ng LH 
R 1—R Se 
Thus f S 
Re Be Rs 1— Ro | e If R? | (92) 
R, R? I—R I—R aa r 


Employing Equation (64) it is now possi sermi : 
S possible to determine the act . 
at any value of the (time dependent) ratio R. er 


ad? — 2 Bar a. | R "1 R? 
R, RS 1—R ie | 


RTL) 
a? 


ee ee 
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From considering Equations (63), (64) and (65) it can be shown that 
| ar 1 el 
=, In rer]: atte, Dee He Ma, re ( 9). 


In the above equation, a’, e& gark and k have already been determined and &, 15 thus a 
function of r alone and can be found for any point in the cylinder wall. 
The maximum value for &, Will always be &,, attheinner radius. Rearranging Equation (92) 


gives 
ea R, 1— Rs Re 
el le er 


0.5 1.0 . 
TIME PARAMETER L>t Lt,=1.33 


Fig. 5. Strains at Bore versus Time for R, = 0.25andn =1 


Since e, and k are now known, the radial strain is, from Equation (1) 
FE) U Re a EINS 
Values of &,4 &a, and k determined from the above equations for the case of n=1 and 
R, = 1/4 are tabulated in Table I and plotted in Figure 5 against L - t. Table I illustrates quanti- 


tatively the extremely rapid increase of strain at times close to the theoretical failure time. 
The strain rates may also be determined in terms of the ratio R (and hence as function 


of time). 


From considering Equations (82) and (87) the longitudinal strain rate becomes 
Rs Ti Ri ni R Ki eE 
Ben 10 allen SL 
From Equations (60) and (80) 
BEDAET: 
I ER Dr EN) 
(&2 kayu2 IR (99) 


at)" u NER TOD. 


ri 


The maximum value of & will always be at the inner wall and 


f 12 11/2 
1 ER nn ION 
u 
Similarly, Equation (36) can be rearranged to give the tangential strain rate 
ee nain.....0..... dom 
? 22 \ FR/J 


33* 


496 F. P. J. Rımkorr and J. R. Luk, Large Strain Creep of Rotating Cylinders 


.— 


Equation (15) gives the radial strain rate 
a ey ee Be “oc 10E 
The creep strain rates at {= 0 can be found by observing that R= R„«=aandb’=b. 
Equation (98) then becomes 
ee. Be N ee a A (103). 
z 1—R 
Equation (100) becomes 
. [a 12 ) in 2 (104) 
PRSAEBBER 7% ad, BEER EN SE E 2 
y le r Si 
Equation (102) becomes : 
a 


SR ) | u a 
| 5 alt! 
Following the pattern as established, the stress distributions over the cylinder wall will 


also be expressed as functions of R. 
Equation (42) may be written as 


_— kim yalm dy 08 


In the case at hand, n = 1 and integration gives 


—E > er (107). 


Be — 
"  y3B 
Now 


After appropriate substitution of Equations (108), (80) and (98), Equation (107) becomes 


R; 77 1—R “| R | 1 | & 0 w*? AR e 
bl 1-1 Ri], a Pen rer 
Remembering that whenn =1 
o 1 
nn ae Le en, (110) 


Equation (34) can be modified to give 


L Il ri R Y 1 | = Bd 
— u 3 Q 
. il R, TzerRt TR I+ m 5 ( a”) Alt 
From Equation (33) it can be shown that 
0, = | L |l— oe R N 0 w* ‚2 2) 
’SWBTER, a 5 N a 


The significant true stress o can be found by employing Equation (8) but if & has been 
determined across the wall, an easier method is to use Equation (22) in the form 


€ 


ee eh En A 


For an analysis of creep stresses at {= 0, Equation (21) can be shown to be 


0 w* a2 b? 
0.,= 5 @ + 2) m B= r) 3 E (114). 
Equation (111) becomes 
KO) am 
- : @ ET E er r) Hilo: 
Equation (112) becomes | 
o we: [/b? + a2 : 
are 77) TE re Bee hie 
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Equations (114), (115) and (116) will give the initial creep stresses in a cylinder made of 
any material without the restriction that n be 1. While the initial stresses may be the same for 
any material, the creep rates will not be as they are determined by the values of B and n in 
Equation (7) BB... RE 


The distributions of strain rates, stresses, and strains across the cylinder wall were computed 
for an example and are shown in Figures 6, 7, 8, 9 and 10. In the example the following data 
were used: B = 18 . 1012 (in21b"h),n =1, 0 © = 6000 (Ib in”), a=1(in,b =4 (in). 


18 10% 100 000 


16 


TRUE STRESS (1b/in?) 


CREEP STRAIN RATE 
fe) 
o© 
[o} 


Fig. 6. Initial Distribution of Creep Strain Rates across the Cylinder Wall Fig. 7. Initial Distribution of Stresses across Cylinder Wall 


ED 
[:} 
Kal 
Fin 
a 
10 - 10°0_ S | 
a 
4 5 
8 E 
ES ii 
3 50.000 
>} 
[>} 
a 
11 
„+ 
8 
[>J 
[77 
5 (in) 
0) ER 
-10 - 10° 


Fig. 8. Distribution of Strain Rates across Oylinder Wall at 35000 hrs Fig. 9. Distribution of Stresses across the Cylinder Wall at 35000 hrs 
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7. Solution forn +1 
Equation (77) stated that 


5 en u u a re ea ne 


As mentioned previously, it is necessary to develop a relationship between either ya? en 
or Yß?— 1 and the instantaneous radii ratio R, so that Equation (77) can be integrated. 
{ It is possible to eliminate one of the 
two integrals in Equation (62) for it can 
be shown that integration by parts of the 
latter integral and appropriate rearrange- 
ment will produce the following inter- 


0.7 


0.6 


0.5 ; dependency relationship. 

0.8 ß „ 

= i "a 1 ji dx " we; 
Mr: zwi 
0.1 5 (in) Fz 


Fig. 10. Strains across the Cylinder Wall at 35000 hrs 


Substituting the above into Equation (62), rearranging and simplifying gives 


B 
1 3 AIME 27 er 
v4 ___ 2yönle eier) RER 
a1 era —-wP=l)+ee ie 
When integrated by means of series, Equation (118) becomes 
bi 
1 n a Br = R. 1 24 
sen) _ aysle "yet m ng 
& @WAPa—2an) Byaya—ıyR-i)+nVa-ırya) 


Since x and f are always negative, it is advisable for simplified computation to rearrange 
Equation (119) to 


& 


1/n = 
B ayR®— i) 
= = — — — _ ; 120 . 
eye tyRld+nVe—TH RN. 
If Equation (118) cannot be integrated directly, Equation (120) can be employed in the 


following manner to relate Ve — 1, yp? -1landR. 
From Equation (74) 


«a |n 


N (a ne x.) 2 /3|# Te 
— 2A) —2An) 


YR—i= RyyaT yes Na 


By investigating R over the range O< R< 1 the entire range of possible cylinder wall 
thicknesses during creep will be explored. 


For a particular value of R, a can be estimated and / ß? — 1 determined from Equation (121). 
These values can now be entered into Equation (120) and if it is satisfied, the correet values of 
Ya — 1 and Yß?—1 will have been found. By following this procedure for other values of R, 
values of Ya? — 1 and YR—1,&and P can be determined and plotted in a manner similar to that 


1 u nd 3 u a EZ TEE 
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illustrated for the case of n = 1 (Figure 3). There would be a plot similar to Figure 3 for each 


value chosen for n and by appropriate curve fitting techniques, Ya® — 1, YB2 — 
be found as functions of R. Ä ‚ ae a ii 


From Equations (75) and (77) 


2dR 


ee er er En ; 
YERA— Ryya—ı (122). 


Knowing yo? — lasa function of R makes it possible to integrate Equation (122) to find k 
and e* as functions of Rand R,as was doneforn=1. . 


Returning then to Equation (85) and rearranging it, gives 


ß 
Rn EN 
ek ja In = Bm) ....... (2). 


[3 


Solving for £ by integration gives 


I s RT 8 a n 

ur eo. n N 

E eu 2 2 Ver EdiE er 
3 BUr 0 02 ( — a?) Ye—]1 


R, 
Employing Equation (118), the above equation can be expressed as 
OR 


| —4n fl Pr Ya—I—or YR_1 


n 


enk dk (125). 


ea) ) laysye—ıyP—1)+n(Y@—ı+yR— 1) 


Ro 


Since everything inside the integral sign has already been determined as a function of R, 
substitution of the various quantities, simplification and integration will give the instantaneous 
radii ratio R at any time /. With R known, the analysis can then be accomplished by following 
the steps outlined for the case n = 1. The creep failure time f, is obtained by setting R= lin 
Equation (125) and is 


f -5| = 
NET: 0 w2 (b? — a?) 


ayaye 1 yRZT Hrn ya ymn)| © ak (126) 


an [/ | Bm+Un vor = 1 + aatı)m Yp? = 1 n 
RK 


8. Conelusion 


The problem of conducting a creep analysis for a rotating cylinder made of a material 
creeping according to an n-power relationship between strain rate and stress has been studied 
and the complete solution for the creep failure time, deformations, strains, strain rates and 
stresses has been arrived at for the case when the creep exponent n = 1. The theory developed 
is perfectly general and by following procedures which have been outlined, it will be possible to 
obtain similar solutions for any value of the creep exponent and hence permit the analysis of 
eylinders constructed from any material exhibiting creep. 

The theory which has been developed is based upon the best available information at the 
present time regarding creep and the behaviour of multiaxial stress systems. The true stress has 
been used and the use of the true strain results in the applicability of the analysis to strains of 
any magnitude from zero to infinity. To illustrate this particular point, a cylinder was chosen 
for analysis which underwent a change in the bore radius from 1.0 inch to 1.86 inches at the time 
investigated. 

It should also be mentioned that the theory was developed in such a manner as to cover 
the entire range of possible cylinder wall thicknesses and that no simplifications were employed 
to treat a particular wall thickness range. The question of whether the cylinder under analysis 
is thin-, thick-, or very thick-walled thus becomes of no consequence as the solution is exact for 
each case. There are also no restrictions as to the magnitude of the inner and outer radii, the 
material density or the speed, and when the theory is extended to all values of n, there will be 
no restrictions upon the material itself as long as it is isotropic and homogeneous and as long as 
its creep behaviour can adequately be represented by an n-power relation. 
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Influence of Large Amplitudes on Flexural Vibrations 
of Elastic Plates 


By NoBoRU YAMARI 


Für die nichtlinearen Biegeschwingungen dünner rechteckiger oder kreisförmiger Platten werden Nä- 
herungslösungen angegeben, wobei verschiedene Fälle von Randbedingungen betrachtet werden. Dabei wird 
der Einfluß großer Amplituden auf freie und erzwungene Schwingungen klargestellt. 


Approximate ‚solutions for the nonlinear bending vibrations of thin plates are presented for the cases of 
rectangular and circular plates subjected to various boundary conditions, and the effects of large amplitudes 
on both the free and forced vibrations are clarified. 


Haxonnmrca IPHÖNMKeHHOe pemieHnHe HEJIUHeÜHLIX BHÖPalmaÄa TOHKUX HPAMOYTOJIBHBIX U 
KPYTOBbIX IJIACTUHOR C PA3HbIMNH KpaeBbIMH ycıoBuamu, BpiacHnerea IhherT 60NBIHMX am- 
IIIUTYA MIA CBOÖONHBIX U BbIHY+KIIEHHBIX Kolebannü, 


1. Introduetion 


A variety of important problems of structural strength and stability of plates, arising in 
modern aircraft construction, can not be adequately analyzed on the basis of the classical theory 
since the plate deflections experienced are not small in comparison with the plate thickness. 
In this case, the interaction between the membrane stresses and the curvature of the plate must 
be considered, which leads to nonlinear terms in the equations of equilibrium of the plate element. 
Among the theories taking into account the influence of large deflections, the most widely known 
is the static theory of von KArMAn, based on which various problems on the large deflection and 
postbuckling behavior of thin plates have been successfully investigated. 

On the other hand, the dynamic behavior of plates subjected to intense acoustic loading 
has become a problem of increasing importance since the advent of turbojet engines, and several 
researches, both theoretical and experimental, have been conducted. For example, HERRMANN [1]!) 
proposed the nonlinear plate theory of motion corresponding to the dynamic analogue of the 
von KArmäÄn theory, while ERINGEN [2] made a comprehensive study of the vibration of vis- 
coelastic plates taking into account shear deformations as well as finite deflections. Recently, 
a new theory which is more general than that of HeRRMANN was developed by TADJBAKHSH and 
SAIBEL [3]. Although the basic equations have been thus established, their applications to 
practical problems seem to be confined to only the special cases owing to the difficulty of the 
mathematical treatment. That is, Cav and HERRMANN [4] studied the influence of large ampli- 
tudes on free flexural vibrations of rectangular plates with hinged, immovable edges, by applying 
HERRMANN’s theory. Approximate solutions for the nonlinear response of rectangular plates to 
sinusoidal acoustic pressure were obtained .by KIRCHMAN and GREENSPON [5] by using the static 
load-deflection relations previously obtained. Experimental researches on the nonlinear forced 
vibrations of rectangular plates and rotating circular discs were made by LassıteEr et al [6, 7] 
and Tog1as [8], respectively. 

In this paper, the extension of the work of Cuu and HERRMANN is attempted. That is, 
based on the simplified plate theory of motion corresponding to the first order approximation 
to that of HERRMANN [1], free and forced vibrations of both the rectangular and circular plates 
are investigated under various boundary conditions. The corresponding static cases have been 
also treated for checking the accuracy of the solution. It is ascertained that the results here 
obtained compare favourably with those previously obtained. 

In what follows, the rectangular and circular plates are treated separately. In each case, 
the basic equations and boundary conditions are first described and then the method of solution 
is given. Finally numerical solutions are presented and discussed. 


2. Nonlinear Vibrations of Reetangular Plates 
2.1 Basic Equations and Boundary Conditions 
Under the assumption that the effect of both the longitudinal and rotatory inertia forces 
can be neglected, the basic equations governing the nonlinear vibration of plates subjected to 
sinusoidally varying pressure can be reduced to the following set of equations?) 


V:F=E(Wy — We. W.y,) EEE uhr 
Lw, F=DV:w+ohw„—h(FyyWas + FıezWwyy— 2 FayWa,) —pceosot=0 (2). 


1) Numbers in brackets refer to the references at the end of the paper. 
2) Subscripts following a comma stand for differentiation; e. g., 
way = &wjox day, Fyy = ®Fjöy’ ete. 
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i nd to the first order approximation to those of HERRMANN [1] as stated 
en equations, x, y are Cartesian coordinates; ıw denotes the plate <a 
tion, h the plate thickness, o the density, £ the time, p cos ot the exciting pressure, and ihr 
E h®]12 (1 — »?) the bending stiffness where E and » are Young’s ra ans EISER EM 
respectively. The operator 7? is the iterated LAPLAack operator 0*/dx* + 2 9*/0x? dy? + 0*Joy*. 
Further, F is the Amy’s stress function defined by 


BPy=Gs Far-tp Po U a . (8), 


are membrane stresses. 
RE the boundary conditions for the problems here treated, Fig. 1. There 
are boundary conditions for both the lateral displacement w and the stress function F. We will 
designate the former as supporting conditions and the latter as 
g stress conditions. 

For the supporting conditions, we will consider the follow- 
ing two cases, i. e., I: All edges simply supported, and II: all 

edges clamped. Thus we have _ 


I:v=w.+r/0,=0, v=Wy+tr’D.=0,t (4). 
B| I:vw=w.=0, v=w,=0 
Concerning the stress conditions, we will treat the following three 
q cases, i. e., (a): All edges stress free, (b): all edges immovable 


Fig. 1. Dimensions and coordinates of tho and (c): all edges movable. These conditions may be expressed 
rectangular plate as follows: 


= +aj2, y= + b/2, 
(a) : N F3z=Faoy-0  # ADaE 
(b) : nr, 0, D=F,-0 RE) 
(oe Doz 0. u consin P,=F.y=0,;, v= const. a 


In these equations, u and v are mid-plane displacements, while P, and P, are resultant forces, 
in the x and y directions, respectively. These can be expressed by the following formulas. 


1 il 1 1 
u = Eu dc v -/[ (F 2. —r Fu) 5 w’,|dy (6), 
0 0 


vj2 a/2 
Ps=hujofyydy: Pak) Pascal Er 
—bj2 —a/2 
As seen from Equation (Sc), we mean by “movable edges“ the edges which are kept straight 
by a distribution of normal stresses, the resultant of which is zero. 
The problem consists in determining the functions w and F which satisfy the basic equations, 
together with the prescribed boundary conditions. 


2.2 Method of Solution 


For obtaining approximate solutions, we first assume the deflections in Cases I and II as 
be w = ht) cos (n/a) x cos (r/b) y, a 
11% w = hit) cos? (r/a) x cos? (r/b) y 8). 
where y(t) is a function of time alone, which is to be determined later. These expressions obviously 
satisfy the related supporting conditions. Here it will be noted that the maximum value of Ye) 
coincides with the maximum defleetion w„ divided by plate thickness h, which will be denoted 
by wm. Thus 


(je Fl nn TE Den en (9). 
Substituting Expressions (8) into Equation (1) and letting 
a Re ei 
the results may be expressed as 
ne BE Dee 2 R 
ET 22 Dog Non Yan» PN 2, an ze er 


where we have used the abbreviations 


X. =co(2pzja)r, Y.=es@gnlb)y.........(2). 
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The general solution, f, of Equation (11) is given by a sum of the complementary function, fs; 
and a particular integral, f,, i. e., 


er EL) 
and /, may be expressed as 
h=ERZ Zoo Yan: RL 
where 
va = P? Dy./16 (p? + P? 0 B=ajb ee en Sr (19) 


Performing calculations, the nonzero coefficients ,, in each case are obtained as 
I: 91 = —1/32P, 90= — B?/32, 
1: 91 =—-132P, 9-=—P32, Pr = — BR/16 (1 + 89%, 902 = — 1/512 9°, * (16). 
= P512, Pe = PHP, 9 = — BR]32 (4 + P2) 
It will be easily seen that /, thus obtained is an even function in x and y, with vanishing shearing 
stresses along the boundaries. 
Noting that /, must have the same properties as stated above, we express /, as [9] 


1 1 
a 


A, 
- Eh2 — = (sinh je “g cosh 7) Ai 
n= E (sin 
inh Z" cosh RT 
B ß Ri 
9n 
ante = “ cos nr 5 
ß a 
Bı | 3 n 
+ 
n? B? (sinhna ßcoshnzß +nnß) b 
2n: r 
En DE zsinhnaßsinh Tr os a a a AR el, 


where p,, Py, A, and B, are unknown constants. Substituting the expressions of w and F thus 
obtained in Equations (5b) and (Sc), we obtain for Case (b) 


1+Pß ®ER ®+v ®ER 
A DB, =0, PD: >= (O8. Ten C a2 N) Py — (5 BT 2 (18), 
where values of C are 1/8 and 3/32 for Cases I and II, respectively, while for Case (c), we have 
ao dc pee DER We 17 re Be 


For Case (a), we first observe that 


n 


An®Eh2 
Vo.yule=+a/2 = Pz Eee rar Ir. A „ed 


a? = 
472 E h2 
a2 


2 rn. 21.420), 


[Ar Xn + 2" DB, &(nB, = 


Vo... lu= +072 — Pu = 


SM 


where 
(sinn yn —y mcoshyn) coshyaz + yarzsinhyasinh yrz 


ei . (ab). 
v2) sinn yrrcoshyan + ya 2 
Expanding the function &(y, z) in a FOURTER series in the interval —1<2z= I, we obtain 
3 4(— 1)” y m? sinh? y 
c zw z h == I = 22 ’ 
m ehe ren Y* + m) sinhymcoshyan ty er 
and then Equations (20) become 
4m®Elh? 2 h2 ae \ er 
Uolensun = De EB + 2 Nr Annlgen]| Yan. | 
N m . (23). 
4n®E h2 : 
Uoaeclu= zur = Pu, 2 [Ant 2 1)" Bmnlm Br] Kon 
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Therefore, in this case, we have from Equations (5a) the conditions p,; = p, = 0, together with 
a set of the following equations for the determination of A, and B,: 
Ei 
Ba+ EN Annlgen)-— MR Nr Een 
m 
An + & 1)" Bmn(mß, n) = —n? - 1 Ina 
m 
m,nal,233..' 127% 
To facilitate the caleulation hereafter, we rewrite the function f, for Case (a) in the following 
double FOURIER series a 
ei El Z WaXz» Yag» 
p 


p(— 1) e, sinh? 22 


aß a a. sa ec Benin in! (3), 
ee To ER: DE Dar F sinhgmBeoshqußtguß "| - 
a an Sen 


pqg=0,12... 


where, = 1/2andg =& =: = 1. Hence, the stress function f in each case may be generally 
expressed as 


1 1 h 
I=-FWPO=ZPWP tg +ERZE PpaXap Ya... - > (26). 


It is evident that for Case (a), pz = Py = 0, Dya = Prag + Pp4; for Case (b), p, and p, are given 
by Equations (18), ®,, = %».; While for Case (ce), p, = p,y = 0, and ®,, = Prr- 

We have thus obtained the expressions for w and F satisfying the boundary conditions as 
well as the first of basic equations. With these expressions, however, it is obvious that the second 
of the basic equations, Equation (2), can not be exactly satisfied. Hence, instead of Equation (2), 
we impose the following requirement on the functions wand F 


a2 b/2 
JS HwFwdidy=0.....ceeeeee: (27), 


above condition implies that the virtual work, as a whole, vanishes at every instant. It is easy 
to see that, for the static case where y(?) is an unknown constant, this method reduces to the 
one-term approximation of the GALERKIN’s method. 


Inserting the expressions for w and F in Equation (27) and integrating, we obtain the 
following nonlinear differential equations in each case. 


{ 
from which we will obtain an ordinary differential equation for the determination of ylf). The | 
| 
| 


n:Dh 2 h2 

I: + + Ay HT + Rp) 

2 rn“ E h4 16 

En Pur d)y=zpesot, 
28). 
IT: R2 16r!Dh > Bra hr 4 h? = 
° ee yın er +2R%+3P)y+ Ew) (Ps + B? p,) y? 
32 4 E h4 


1 1 16 
a (0. + Po+ Dı+ 5 Dart 5 Dar + Dat Duo) y’ =„pcosot 


PT RB RER equations for y(l), both the functions w and F will be determined and the problem 
is solved. 


2.3 Numerical Solutions for the Square Plate 
For the square plate, i. e., 8 = 1, Equations (28) become 


j J 4 E h% 1 16 
I: ya a samt] ponot, 


“ 


1% h2 
opt zu mrtE 


| 2 


16 
ayer\ w| = post 


m. K1, Influence of Large Amplitudes on Flexural Vibrations of Elastic Plates 
In these equations, values of £, are 0.06492, (3 —»)/8 (1 + v) and 1/8, while those of £, 
are 0.14903, 0.16656 + (0.14063)/(1 — ») and 0,16656, for the stress conditions (a), (b) and (c), 
'  respectively. In calculating A, and B, for Case (a) by Equations (24), we have taken n from u: 
1 to 8. Now we will first treat the static case for checking the solution, and then proceed to 

evaluate the effect of large amplitudes on the free and forced vibrations. 


(i) Static Deflection 


Letting Yu = 0 and cos = 1 in Equations (29) and recalling Relation (9), we have the 
expressions for the static maximum deflections in each case as 


28 1 
1: ln mn 
E H4 16 3 (13) "m + 51 Wm| » 


s 2 
ib 2 = en, mi 
E H4 2 3(1—»2) Om + Co wm 


Fig. 2(a). Maximum deflection of square plates with Fig. 2(b). Maximum deflection of square plates with 
simply supported edges; v = 0.316 clamped edges; » = 0.316, not otherwise specified 


The numerical results are illustrated in Figs. 2(a) and 2(b). In these figures, (a), (b) and (c) 
specify the stress conditions previously defined, while (I) designates the corresponding linear 
solutions [10]. Some of the problems here considered have been solved more accurately by Levy 
[11,12] and Way [13], which are shown together by dashed curves. Poısson’s ratio » is taken 
as 0.316 except for the WAy’s solution where » = 0.3. It will be seen that in spite of a rough 
approximation, the present solutions are in reasonable agreement with those previoulsy obtained. 


(ii) Free Vibration 
The equation for the free vibration of square plates is obtained in the following form by 
putting p = 0 in Equations (29): 
Vrr = 2 y + u 1° = 0 PIERRE (31). 


For vibrations with small amplitudes, the term u? y? can be neglected in comparison with 
2» and the equation becomes linear. The corresponding natural circular frequency is given by o, 


which becomes (2 r2/a®)Y D/o h for the simply supported square plate. This value precisely coin- 
cides with that from the linear theory [14]. For the clamped square plate, x becomes (37.22/a?) 


YD/o h, which gives 3.5%, larger value than the accurate one [15]. 


Assuming that the amplitude of the vibration y is equal A, and it oceurs at{=(), the 
solution of Equation (31) will be given in the form of an elliptic cosine en 


Bene. ro. EEE EAN, 


h 
Se A= (0? + u2 A2)I2, BEAT u en ra 
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The period T of en(At, k) is 


TZAÄARM 2 u 20 ne a ee 
where K(k) is the complete elliptie integral of the first kind. The corresponding linear period, 
Ir, is ? 

- RT... (35), 


which may be deduced from Equation (34) by observing that as A —>0, A —o and K-npß2. 
The ratio T/T, is given by 
e3 = 2uaR ” . . . . . . * ” . . . u . (36). 
T, yo? + u? A® i 
In Figs. 3(a) and 3(b), this ratio is plotted versus the dimensionless amplitude pn assuming 
that » = 0.3. The meaning of the symbols in the figures is the same as before. It will be seen 


that in each case, the period decreases with inereasing amplitude and that this tendeney is most 
pronounced for the plate with simply supported, immovable edges. 


06 
U Ba ee 3 Fe de dr 00 UT ZUBI SR er 
Fig. 3(a). Period of vibration of square plates with Fig. 3(b). Period of vibration of square plates with 
simply supported edges; » = 0.3 clamped edges; » = 0.3 


Chu and HERRMANN [4] have studied the influence of large amplitudes on free vibrations 
of rectangular plates with simply supported, immovable edges, by using a perturbation method. 
In this case, we obtain from Equations (28) 

Dh 3n:Dh 
++ E-MU+NH+HAIP=0 N, 


A 


7* 


oRyu+ 


which is ascertained to be precisely equivalent to the fundamental equation, Equation (20) of 
Reference [4]. 


Here, it is to be noted that the following approximate formulas [16] give satisfactory results 
for values of A less than |a/a|: 


Hi ( 3 12 2 
Fa re ( a ö Er a) i 


| 


To N, sr 
3 uUAS EEK: j 
yil)=Acosat+ 35% (cos 3a’ t— cosa’) 
(tii) Forced Vibration 
Equations (29) may be generally expressed as 
Yyarayrı Ay = Psalm. 2 2 Sr 
which is the well-known Durrıne’s equation. The approximate stationary solutions are given 
by [17] 
w* 3 8 r 
a ia ER 
a2 m 4 02 - A a8’ | 
2 AB LEN 


vll) =Acoswt + 


i n rei 
32 |2 Em aA 
2 (x + A 7) 


where A is, as seen in the second of these equations, the amplitude of the component with the 
cireular frequence w. 


ER 
‚ 
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Fig. 4. Frequency response of the square plate with simply 
supported, immovable edges; v = 0.3 


pP 


3. Nonlinear Vibrations of Cireular Plates R; 
3.1 Basic Equations and Boundary Conditions 


Transforming the preceding basic equations, Equations (1) and (2), in the polar coordinates 


and noting the axial symmetricity of the fundamental mode of the vibration of circular plates, 


we obtain 
| E | 
: Wer EN RAS TEN, 
‚ h 
L’w, F)= DViw + ohw.— (FR, We = peswot=0 ...... (4), 
where 
1 en =) 
a 
V r dr dr Er ee ae 


- Both the functions w and F are evidently independent of 0. The connections of the membrane 


stresses with the stress function and midplane displacements are 
given by '& 
1 BE Ei v 
%"=—F,= (u. E Wr zer 
De‘ Les TR 2 r 
(44), de 
E 
n=F rr 


a a 


1 v 
e Ltr tn Wr 


where u, is the radial displacement. Tangential displacement and 
shearing stress are identically zero from the assumption of axial 
symmetry. 2R 

Now we will state the boundary conditions for the problems #ie-5- Dei alnlenegie 
here considered, Fig. 5. For the supporting conditions, we consider 
the following two cases as before, i. e., I: elY supported and II: clamped. Thus we have, at 
TIER: 


8: “ —0 
T- en ala e ba den Bam Du y BirinT (AB) 
IT: Pa a 


Asto the stress conditions, we will treat the two cases where the edges are free from stresses and 
constrained immovable. Designating these cases by (a) and (b) as before and noting, from Equa- 


tions (44), that 
N... en. nu (db), 


2 


E N. Yanakı, In Inf fluenc 


we have the following SSRUHAnNE at es Br 
@: li Ve | 


Ba 
(b) : er Mir or 


Be. The problem consists in finding the solutions w and r of ‚the b 
Re given boundary conditions. , 


x | Kr 
2 3.2 Method of Solution MER E ı« 
j 'To obtain approximate solutions, we assume the deflections in RPRSRABe BL .; 
- 


Fi | w=hy|i+G(E) +2] ROT 2 2.68 


f where y(?) is a function of { alone, while C, and C, in each case are given by 


v2 GC = 
IT: Bess U | 
It is easy to see that the above expressions satisfy the related boundary conditions, and that } 


Bm ern. : (0) 8 
Substituting Expression (48) into Equation (41) and letting 


FE eg EEE 6) 


the solution R the resulting equation may be expressed as j 
4 6 8 
s 2 =) = (2) | ER 
erlag RAR 62), 


where C, is an unknown constant. Inserting this equation into Conditions (47), we obtain (C, 
in each case as 


GG yRBAHICaH2N, 
HEBEN aHE GG +20— ng] 


With these expressions of w and F thus obtained, we impose the following requirement 
instead of Equation (42) as before: 


. (53). 
(b) : Gg=— 


1 Liw Fwrd=0..,.1202 nn EE (54). 
Taking v = 0. 3 and performing integration, we have the differential equations for yff) in each 
case as 
2 Et / F 3 
ie oRyu+ Ri (2.242 9 +, y°) = 1.560 p cos wt, 
x Are (55), 
R Eh > a 5 
E II% oh Yatzr Ey HP) = peoswl 


where values of £, are 0.591 and 4.148, while those of & S, are 1.429 and 4.602, for the stress con- 
ditions (a) and (b), respectively. 


3.3 Numerical Solutions for the Cireular Plate 
(Ü) Static Deflection 


Letting 9... = 0 and cos ®! = 1 in Equations (55) and recalling Relation (50), we have 
the following equations for the maximum defleetions in each case. 
1.438 w.. + 0.3789 io}, for Case I(a), 
p Rt 1.438 wm + 2.659 13, for Case I(b), 
Eh 9.861 wm. + 0.8571 ı$, for Case-Tlla),.f m Se 
9.861 0, + 2.761 1}, for Case Il(b) 


o© 


(=) 


f 
f 
en 
% 


Ai 2 


BURN 


BEN 75. 


LIFE 
arr 


LA 2 3 
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En corresponding approximate solutions obtained by TımosHenko [18] for the same value 
of v are 

1.437 w,. + 0.3764 ı, for Case I(a), 

1.437 wm + 2.661 wm for Case I(b), 


Eht 15.8480, + 0.853880, for Casella),f 9 


5.848 w,, + 2.754 w), for Case II(b) 
which may be seen to be in excellent agreements with the present solutions. Further, solutions 


of exact nature have been obtained by Sriprzs and Hausrarn [19] for Case I(a), while by WAY [20] 
for Case II(b), which are ascertained to be in reasonable agreement with the present ones. 


en) 
2 


Fig. 6(a). Maximum deflection of ceircular plates Fig. 6(b). Maximum deflection of eircular plates with elamped 
with simply supported edges; v = 0.3 edges; v = 0.3 


Numerical results are illustrated in Figs. 6(a) and (6b), together with the previous solutions?) 
as well as linear ones [21]. 


(ii) Free Vibration - 
Putting p = 0 in Equations (55) gives the equations for the free vibration having the same 
form as Equation (31). The linear natural frequencies, i. e., &, in each case become 


1.497 ıB £: for Case I, 
a 
6, EE (58), 
E e 
I5.125 Ir VE for Case II 
which are found to be 0.2% and 1% larger than 


the corresponding exact one [22], respectively. 

The influence of large amplitudes on the 
period of free vibration can be obtained by the 
same procedure as before. The results are pre- 


sented in Fig. 7. It may be seen that this effect 02 
is greatest for Case I(b) as in the preceding case 0 04, 3085742 1,6 2 
for square plates. Fig. 7. Period of vibration of eircular plates; » = 0.3 


(iii) Forced Vibration 
Equations (55) are obviously in the same form as Equation (39); hence, the nonlinear 
response to the varying pressure will be obtained in the same manner as before. 


Acknowledgment 


The author wishes to express his thanks to Professor Emeritus S. HısucHı and Professor 
F. Numacuı, Tohoku University, for their invaluable advice during the course of the present study. 


3) The results obtained by TımosHENKo and Srıpres et al are not shown in the figures, as they are not 


distinguishable from the present solutions. 
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KLEINE MITTEILUNGEN 


Die Vergrößerung des Halbmessers eines 
durch Pole begrenzten Konvergenzkreises 
von Potenzreihen 


Untersucht werden hier nur Potenzreihen, deren 
Konvergenzkreis durch endlich viele Pole begrenzt 
ist. In der Praxis kommt es oft vor, daß der Wert 
einer durch ihre Potenzreihe ausgedrückten Funktion 
in der Nähe der Pole, wo die Konvergenz stark ab- 
nimmt beziehungsweise ganz aufhört, zu berechnen 
ist.. Zweck der vorliegenden Abhandlung ist es, den 
Radius r des Konvergenzkreises zu vergrößern. 

Zu diesem Zweck sei folgendes Verfahren emp- 
fohlen: E 

Von der Funktion wird die Summe der zu den er- 
wähnten Polen gehörenden Hauptteile [1] abgezogen- 
Diese neue Funktion besitzt auf der Kreislinie mit 
dem Radius r keine Singularität mehr, naturgemäß 
auch nicht innerhalb des Kreises. Infolgedessen ist 
der Konvergenzradius R der Potenzreihe der neuen 
Funktion größer als r. 

Die Funktion wird mit f(z) bezeichnet. Der Mittel- 
punkt ihres Konvergenzkreises sei c und sein Halb- 
messer r. 


Ihre konvergierende Potenzreihe ist also 
fe) = u +1 @—0)+..-+m@—c)" +... (1). 
%®—e <r 


Die Pole sollen sich auf den Punkten »,s,...,t 
der Kreislinie mit dem Radius r befinden. Ihre Ord- 
nungszahlen sollen 7,0,...,r sein. Die Haupt- 
teile der Funktion sind: 


P, P, 2 


u ee een, a per 
A S, 8, N 
Ys(2) = 2 8 — (z — s)? IT 2% Sk (2 — s)° ’ 

er; Ta Dosen 4 
9:2) tree ak 


Die Bezeichnung der Summe der Hauptiteile ist 
plz) = Ppl@) + 92) + + il). 


Diese Funktion ist rational. Singularitäten hat sie 
nur in den Punkten p,s...t der Kreislinie vom 
Radius r, sie kann also um das Zentrum c in eine 
Potenzreihe entwickelt werden: 


th) + + Bett... 
Belr<r: 
Diese Gleichung wird von Gleichung (1) subtrahiert. 


f2) — pl) = (ao — do) + (1 — bi) R—e)+: 
+ (an —bn) 2 — ce)" +: 


Wie oben bereits erwähnt wurde, hat diese Funk- 
tion keine Singularität auf der Kreislinie mit dem 
Radius r, ihr Konvergenzradius ist R > r. Somit ist 


f(z) = p(2) + (ao — 5) + (aı — di) ®—c)+- 
+ (en in) (2 —e)" +: R—cl<R (2). 


Die Funktion f(z) wird also als die Summe einer 
rationalen gebrochenen Funktion und einer Potenz- 
reihe ausgedrückt. Sie ist in jedem inneren Punkt 
des Kreisgebietes mit dem Zentrum c und dem Halb- 
messer R unter Ausschluß der Polpunkte 9,s...t 
konvergent. 


Beispiele 3 


1.) Die um den Ursprung (ce = 0) entwickelte 
Potenzreihe der Funktion f(z) = sec z ist: 


1 5 T 

— —z2 erg a, De 
sec 2 ten kI<5 (3). 
Die dem Ursprung am nächsten liegenden Pole der 
Funktion seez sind » => s= 7: also ist r -7 E 


Die Summe der zu diesen Polen gehörenden Haupt- 
teile ist 
1 1 7 
a, = BER 
F] IE, 


Diese rationale gebrochene Funktion kann als 


geometrische Reihe folgendermaßen ausgedrückt 
werden: 
TC 4 16 64 7 
er a8 22 ai, 24 zI a 
= er ® Br, D 


Durch Subtraktion von der Gleichung (3) ergibt sich 


sec2z — Sr Hi ir ar 
4 
5 64 = 


Der Konvergenzradius ist auf R = & rı gewachsen, 
da sich die jenseits r = - gelegene nächste ER 
tät a Funktion see z an den Punkten j” und 
Fragrt befindet. 


Für manche praktischen Zwecke geben die ersten 
drei Glieder genügend genaue Resultate 


3,1416 2.097394 0,0162 (8). 


seo 2 ya 


l<za 


In Polnähe ist der Nenner mit größerer Genauigkeit 
in Rechnung zu stellen. 


Die Vergrößerung des Konvergenzhalbmessers auf 
5n/2 kann auf ähnliche Weise erfolgen und führt zu 
dem Ausdruck 


1 3n 
secz = _— 
u2 - 972 
4 


4 
lan = 
q i BIER 


1 16 16\ „ 
+ +57 at 

5 64 64 
+ 


5 
A re nn 
PER"; ») = M<y® 


t 
2.) Aus der Potenzreihe von =: ‚ ferner aus der 


Summe o(2) = —- - der zu den Polen + z ge- 


AERRFER 
4 
hörenden Hauptteile folgt — ähnlich dem Beispiel 1.) 
tg2 _ 2 RR: 32\ , 
2 Pu s ( ==) a nt“ N 
PT z 
may, 3 
Fr Ar; ”) + | < Dr (6), 
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beziehungsweise 
tg 2 2 
le ‚18943 + 0,00482 2? (7). 
ET Sees rn 
< a 
| gr 


Die nach den dreigliedrigen Formeln (5) und (6) 
berechneten Funktionswerte werden in der folgenden 
Tabelle mit den genauen Werten verglichen, wobei 
die vorigen in Klammern stehen. Die Werte von 
(thx) wurden aus der Formel (6) auf Grund der Be- 
tgix thx 


ziehung —; berechnet. 
UL x 

ige) | 182 Si 
z (sec 2) Bec z (#2) | Fi = = = 
0 1 il 1 1 DE 1 
0,2  1,02034| 1,02034 | 1,01355 | 1,01355 | 0,2 | 0,98688 | 0,98690 
0,4 | 1,08573 | 1,08570 | 1,05698 | 1,05698 | 0,4 0,94987 | 0,94987 
0,6 1,2117 | 1,2116 | 1,14020 , 1,14023 | 0,6 |0,89506| 0,89508 
0,8 1,4357 | 1,4353 1,28697 | 1,28705 | 0,8 |0,82997| 0,83005 
1,0 | 1,8517 | 1,8508 | 1,5572 | 1,5574 |1,0 0,7614 | 0,7616 
1,2 | 2,702 | 2,760 | 2,1426 | 2,1435 |1,2 0,6943 | 0,6947 
1,4 | 5,886 | 5,888 | 4,1405 | 4,1413 |1,4 0,6317 | 0,6324 
1,6 -34,241 |-34,247 |-21,3968 |-21,8953 |1,6 0,5749 | 0,5760 
Anwendungen 


Die kleinste kritische Kraft des Stabes in Bild 1 
hat den Wert 


REEL 


Veh —Z 4 a7, . 
AERO. 
22 a 
EI=const. ei 
Bild1 


Auf den Stab wirkt eine Druckkraft P<P, ein. 
Es wird bei verschiedenen Querbelastungen das Ver- 


hältnis —- und = des größten Biegungsmomentes 


0 0 
und der größten Durchbiegung gesucht. Darin ist 
der Zähler das Moment beziehungsweise die Durch- 
biegung unter Mitwirkung einer Belastung P und 
der Nenner dasselbe für den Fall P=0. 


PR 
Mit der Bezeichnung 2 = va folgt 


Be Pr re 
u 422’ 
daraus folgt 
: n® 
en Eee): 
und 
EI 4 
ae (10), 


Das Verhalten des Druckstabes bei Querbelastung 
läßt sich durch die Funktionen sec 2 und we dar- 
2 


stellen. Deshalb soll als Vorbereitung die oben an- 
gegebene Formel (9) für 2? in die Gleichungen (4) und 
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(7) eingesetzt werden. 


T 
tg 2 2 8\ (1 3% m 
z; m 7% + )+(5 „je 
4 4 
2128| m 
+7 »)i6r 
8 % 
ze EL 3 | “ 8\1 
rl) 
0,81 0,012 , 0,001 _»—0,19 


an ee en 2). 
Fa Tu v + WORIZERT (12) 

I.) Der Stab in Bild 2 wird außer durch die Druck- 
kraft noch .durch eine zusätzliche gleichmäßig ver- 
teilte Last q belastet. 


2 


Bild 2 


Das Biegungsmoment in der Stabmitte ist ohne 


Druckkraft M, = es ql2, beim Vorhandensein einer 
Druckkraft [2] 
„ see z — 1 

22 i 


M=gqgl 


Das Verhältnis der beiden Momente ist 
M 2 secz—1 


M, AN: 
Nach Gleichung (11) ist 
e 
n RR | St 4\ı 
gegeh + = +) 
nach Gleichung (9) 
I: 
amt 
Das Produkt beträgt 
gees—1_32 » |, 32, 5m 321 
a matt =): 
v + 0,032 0,004 
era In na FA 


Das zweite und die nachfolgenden Glieder werden 
vernachlässigt und damit wird 


M »+ 0,032 


. 


‚Somit nimmt in diesem Verhältnis das größte 
Biegungsmoment beim Vorhandensein einer Druck- 


kraft noch zu, wenn die Knicksicherheit die Größe » 
besitzt. 


Ist die Axialkraft eine Zugkraft von der Größe H, 
so wird 
2:8 


„- vor 


. Es sollen noch bei derselben Belastung die größten 
Durchbiegungen verglichen werden. 


Ohne Druckkraft: 
bag 


h=anrt 


Beim Vorhandensein einer Druckkraft P: 


gt (zer) =) [2], 


aPp\2 
somit 
12EI/2 
Mr a a (zen). 
Nach der Gleichung (10) ist 
12 EI 48 
3 PR 5m” 
und nach der Gleichung (12) 


2 yerya2/aa 5m2 .32\1 
Be ev 
Das Produkt beträgt 
2a 1, 2 40.004 
e fo 55 v—l ae 1 
3 IL.) Der Stab in Bild 3 wird außer durch die Druck- 
id kraft noch durch eine Querkraft Q in der Stabmitte 
r belastet. 
I Q l # 
pP] 1 or 
Larrg + 
Bild 3 


Die Biegungsmomente in der Stabmitte haben die 


Größe BR 
= reg 
My vr 2 E) M = D} 2 [ ] ® 
Bei Verwendung der Gleichung (12) ergibt sich 


MO tgeT7 90,19 
ee 


Die Durchbiegungen in der Stabmitte sind 


gr zone. 
h=gs57T I-sP = ) a. 


Bei Beachtung der Gleichungen (9) und (11) 


ee 
ss z 


8 

2 2 8 1 
aan | 7 r RE ;) 
a ww —1 12 2,9 

F. DiscHinGeg erzielte bei seiner Untersuchung von 
Druckstäben auf iterativem Wege dasselbe Resul- 
tat [3]. 

Es soll betont werden, daß die einfachen Formeln 
dieses Beispieles auch in der Umgebung der kritischen 
Kraft eine gute Näherung liefern, da der Konvergenz- 
radius der mit den Hauptteilen ergänzten Potenz- 
reihen größer geworden ist. Überdies tendiert der 
relative Fehler gerade in dieser Nachbarschaft gegen 
Null. Somit erscheint die in der Literatur zu findende 
Vorsichtsmaßregel unbegründet, derzufolge in der 
Nachbarschaft von v» =1 die angegebenen Formeln 
mit weiteren Gliedern, die die Potenzen 1/v, 1/v?,... 
enthalten, verschärft werden sollen. 


9 — 0,014 
v—1l 
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So ist zum Beispiel im Falle von v = 1,01 
M _v» + 0,032 
ee] 


Eine einfache Rechnung zeigt, daß die Vernach- 
lässigung nur — 0,005 ist. 


—= 104,2. 
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Steady Temperature in a Composite Circular 
Cylinder 


The method of finite HAnkEL transforms is used 
for determining the steady temperature distribution 
in a finite cireular cylinder consisting of any number 
of homogeneous and isotropie layers. A considerable 
part of what follows is new. 


1. Statement of the Problem 
The cylindrical body 
(ah "VERS 
consists of n homogeneous and isotropie parts 
Izpo=sr, ,_,22S2%, 
Kelle ann 


with the conductivities A,, respectively. We have 
to find the steady temperature field in the body 
provided that its surfaces z = 0, 2 = z, are kept at 
prescribed temperatures f(o), g(0), whereas the curved 
surface 0 = r has steadily the temperature zero. 
The temperature and the heat flux are assumed to be 
continuous in the whole extent of the body in question. 


Mathematically the problem reduces to finding the 
solution ur = ur(0,2), ISk=<n of the equations 


0, 


ur 1 du, | our _ 
20? 0 00 Pu 
Iso, Wer zen el cken 


0488. (1.3), 


u (0, 0) == (0), Un(O, 2n) TE g9(0); 


U;11(0 zk) == UrlO, 27) ’ 


peo<r,, Leben 
a u 
02 02 Ayrı 


"so<r, ve. 1ZzE SR —I 


uln,d)=0, 42, <2<, l12ken (14). 


2. HAnKEu Transforms of (1.1)—(1.4) 
If we apply the familiar transformation [1] 


h(o,) = H[k(o)] = J o h(o) Jo(0; 0) do , 


an] yo) rd 
h(e) = H- [ko] = Ne) Ten 


(2.1), 


where the sum is taken over all the positive roots 0; 
of the equation 


ka) =U ar sum) 


- 
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then the above relations (1.1)—(1.4) change into the On the whole, the solution of (3.2), (3.3) appears 


equations in the form z 
a _g%=0, 4 _1<2<a, ISkszn (23), A, = fe), 
% 1 pr i 
%(0,0) =f(0), UnlQn zu) = 0) (2.4), Bı = Neo) [— Z(g)) fe) + Kon)» 


%.,.09% 2) = Wo, I1Sksn—1 Ayzıl| _ eo || Fo) Ne) 
zn, du | (2.5) Byrı)) No) ||—flo;) Ze) + Fo) (3.7). 
—-—- A —. 2=%, LÄSksn—l1l Eee 


for determining the HAnkEL transforms 


%.(0» 2) = Hluxlo; 2)] 
= fo url, 2) Juloio)do, ISkszn (2.6), 
0 


of the required functions «;(o, 2); of course, ‚No 
9(0;) mean the transforms of f{o) and g(o), respectively. 


3. Solution of (2.3)—(2.5) 
The equations (2.3) give generally 
%rl0; 2) = Ar ch &@ —%—1) + Br.sh 0 — 2%, _ı)> 
Isksn (3.1); 


and with these expressions we find from (2.4), (2.5) 
the linear set 


No) =D, chain + > Shoiln» 
Z(0,) == guD, ch 0; „+ ed, sh o; In 


Substituting these expressions into (3.1) gives us 
the transforms %;,(o;,z) of the required functions 
ur(O, 2). en 


4. Solution of (1.1)—(1.4) 


Having found, by means of (3.1) and (3.7), the 
HANKEL transforms %;(p;, 2), it is quite easy to get 
the solution of our beginning problem. This happens 
by means of the second formula (2.1) which gives 


url0, 2) = . = [Ar eh g,(2 — 2, _ı) 
&; 


BR Br sh 01(2 vi Z_,)] Jo(l0i 0) 


2 2 Jia: rn)’ 
A, =fle), AnchQ;ln+ Bnsho;ln = 90:) (3.2), I<Sk<n (l); 
A; = Arche; k + Brshog; kk, the coefficients A; = Ax{p;), Be = By(o;) are deter- 


1sksn—1 
Brrı — Ar(Ar sh Oi I + Br ch Oi I.) , 
Teken ol 


(3.3) 


of equations for the unknown coefficients Ar, Bj; the 
quantities /; in (3.2), (3.3) are defined by 
ver ee RE 
The best way to solving the above system (3.2), 
(3.3) leads probably over the use of the matrix 
analysis — see [2], [3]. Introdueing the matrices 


mined by (3.7) and the sum is taken over all the 
positive roots 0; of the characteristic equation (2.2). 
5. A Homogeneous Cylinder 

As a first example we shall treat the case 

A=%=:-=4A_1ı-1 

corresponding to a homogeneous isotropie eylinder 
0sosr, 0 S2Sz. 

The matrices Q;(o;) become 


BR een h 
| - = " = -— ‘ 
’ 2x(0;) sh 0; 24, ch 0: 2. a Mg ee. 
k 

Kr = Br)’ lsken, and the formulae (3.7) give 

_||ehgilk, she; . (8.8), Ne) =shaim, ZiQ) = chQizn, 

Ar sho;le, Ancho;lk Ar ch (2 — 2, _1) + Brsh i(@ — 2, _,) 

ISsksn—1l 


enables us to write the equations (3.3) in the form 


Kyrı = Mı Kr, Teken-1T, 


and this yields, with account of the first relation (3.2), 


[Xo:) sh 0:(2n —z) + 5(0,) sh 0,2], 
LsSs#S#: 


Introducing this into (4.1) yields the solution of 
our present case in the form 


2 


= 1 
sh 0; 2 


ee 92) = 5 I az, NO) shosln —2) 
a Sun en y2 (eo) sh 0; 2u 


(8.6), 
Que) = Mr M,_ 1... M,M,, eu 


lIsksn—l1l 


Denoting by rd (r, s= 1,2) the elements of 
Q„_—1(0:) we find from the last of the equations (3.6) 


An = ln !o;) ir q a DB 


n—1 
2 a 9% 
Bn = ed (0:) 2= ZN B 
and combining this with the second condition (3.2) 


gives the value of B, in terms of 0; All the other 
coefficients Ax(o;), Bi(o;) follow from (3.6). 


= Jo(0; 0) 
:) sh 0; 2] —$ i 
+ 90;) sh o; 2] Tloır) 


0<seo<sr, 0<z<z(l); 


the sum is taken over all the positive roots 0; of the 
equation, (2.2), J,, Jı are BESSEL functions of the 


first kind and the HankEL transforms f(o,), 9(0;) are 
defined by the first relation (2.1). 2 9ei 


6. A Two-Layered Cylinder 
The only matrix Qı(0;) is of the form 


ER: cho;2, sh 0; 2, 
Q.(B;) A, sh 0; 2, ! A, ch 0; 2 


and the expressions N (0), Ze) irom (3.7) become 
Ne) = sh +(Aı—1)chg;2,shog;k,, 
Zo)= che; + (A —1) sh 0;2,sho;l,. 


Using again (3.7) we find without diffieulty 


N(0;) [41 cho;2+ B, sh 0:2] = fe;) [sh 0; (22 — z) 
+(A,— 1) sho;l ch o,(2ı —2)] + %o:) sh 02» 
N(o;) [Ag ch o,(@ — 2,) + Bzsh oe — 2,)] 
— /ı fe) sh 0;(2%, — 2) 
+ 9(0;) [sh g;2+ (Ay —1) che; 2, sh g,(@—2ı)]; 


and the solution (4.1) appears in the form 


en f 
10 = Eminem 1 [he »—2) ] 


+(A— 1) sho;l; ch o;lzı —2)] + o;) sh 0,2}, 


_2 J0(0: 0) m 
u,(0, 2) = 5 = N) 0:7) {A, f(o;) sh Q,(2. — 2) 


+ 9(0,) [sh e;2 + (A1 — 1) che; zı sh gil2 — zı)]} 
N(o,) = sh o; 23 +(A,— 1) cho,2,5h0;%» 


=%3—2 


(6.1); 


the sums are taken over all the positive roots of 


(2.2); f(o;), 9(0;) are the HAnKEL transforms of the 
given functions /(o), g(0)- 
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Ein hinreichendes Konvergenzkriterium und 
eine Fehlerabschätzung für die Iteration 


"in Einzelschritten bei linearen Gleichungs- 


systemen 
Zur Lösung eines linearen Gleichungssystems 
ee tl), 


dessen Koeffizientenmatrix A = A; + Ar Summe 
einer linken unteren Dreiecksmatrix A, (mit Elemen- 
ten ajk, 7 > k, und sonst 0) und einer rechten oberen 
Dreiecksmatrix Ar (mit Elementen a;r, j Sk, und 
sonst 0) ist, berechnet man bei der Iteration in Einzel- 
schritten, ausgehend von einem beliebigen Vektor x,, 
eine Folge von Vektoren nach der Vorschrift 


Im+ı — Ariyrı * ÄRm + r Auer hen WE 
(m 0)172) 


In dieser Mitteilung werden dazu ein neues hin- 
reichendes Konvergenzkriterium und eine neue Fehler- 
abschätzung, beide gültig für beliebige Matrixnormen, 
sowie allgemeiner als bisher die Fehlerabschätzung 
von W. Dück für beliebige und die Fehlerabschätzung 
von L. CoLLarz für monotone, absolute Matrixnormen 
hergeleitet. 


1. Bezeichnungen 

Der Raum R, der Vektoren x und der Raum M„ 
der quadratischen Matrizen A sind nach Einführung 
einer Vektor- bzw. Matrixnorm [4] Banachsche 


Räume. 
In M,„ wird eine Halbordnung A < B erklärt durch 


Gi djr reell, O <a Sb (,k=1,...,n). 


Kleine Mitteilungen 


Mit |4| bezeichnen wir die aus den Beträgen der 
Elemente von A gebildete Matrix. n 
Eine Matrixnorm heißt: 


absolut, wenn |||A| || = Mel nee (3), 
monoton, wenn aus A< Bfolgt ||A|| < ||B|| . (4), 


zu einer Vektornorm passend [4], wenn 


4 z|| < |[A]| |}e]] für alle Ae M„undallexeR, (5). 


Das Maximum der Beträge der charakteristischen 
Zahlen von A wird Spektralradius '0(A) genannt. 


Die Spektralnorm ||A|| = + Vol#’ A) ist mono- 


ton [5], aber nicht absolut, wie Gegenbeispiele zeigen. 
Andere gebräuchliche Matrixnormen, wie z.B. 


Ill = Max Z jap, |jAll = Max 2 ja 
5 
sind monoton und absolut. 


2. Das allgemeine Iterationsverfahren 
Das Iterationsverfahren 
"m+1 = Lam + 8 
zur Lösung von 


(m=0,1,2,..) (6) 


IN N (7) 


konvergiert genau dann [5] für beliebiges x, e Rn 
gegen die dann eindeutig existierende Lösung von (7), 
wenn o(L) <1 ist. 
Es existiert dann [4] die Matrix (E— L)-!, und 
aus (6), (7) folgt 
Tan (ED Lam m) (8) 
und daraus (mit zur Vektornorm passender Matrix- 
norm) die F'ehlerabschätzung 
em 4ı — || S |I& — Z) ZI | 1 —am|| (9). 
Ein bekanntes hinreichendes 
Kriterium von J. Weissinger [8] für o(L) <1 ist 
(Eee (10). 
Unter dieser Voraussetzung ||Z|| <1 kann die 
Fehlerschranke in (9) weiter abgeschätzt werden: 


L 
ans al og em ml aD. 


3. Ein neues hinreichendes Konvergenzkri- 
terium und eine neue Fehlerabschätzung 


Da die n-te Potenz der Dreiecksmatrix A, gleich 
der Nullmatrix und deshalb 


(B—Ar)(E+Ar +4, + +40 )=E 


ist, existiert (als Summe einer endlichen geometrischen 
Reihe) die zu # — A, inverse Matrix 


(B—A,)"=E+4A, +47 +4---+47' (12). 


Das Iterationsverfahren in Einzelschritien (2) zur 
Lösung von (1) kann daher als ein Iterationsverfahren 
der Art (6) zur Lösung von (7) mit 


= (EI AT)NARIMEERS 9 (13), 
Ele N le 
aufgefaßt werden. Aus (13), (12) ergibt sich 
ZI = |(& +4, +45 + +47") Anl] 
< (1 + |Aull + Az? + + |Azlle=n |Anl 


ı Azll® ,| 
- ||A falls ||Ar|| # 1 
0 | 1 4zl] Az falls || Anl 
W . |Arl| falls Az|| a 
und aus ||Z]| < O nach (10) das 


Konvergenzkriterium: (15) 
C<L1 ist hinreichend dafür, daß die Iteration in 
Binzelschritten (2) gegen die dann eindeutig existierende 
Lösung von (1) konvergvert, 
und daraus nach (11) die 
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Fehlerabschätzung: 
; 10) 
am+ı — ll SI m ml 10). 
m=0,1,2,...) 


4. Bemerkungen zu den Fehlerabschätzungen 
von L. Collatz und W. Dück 
Die Darstellung (6), (13) ermöglicht die Herleitung 
schon bekannter Fehlerabschätzungen von W. Dück 
[3] bzw. L. Corrarz [1], [2] und zwar allgemeiner als 
bisher für beliebige bzw. für monotone, absolute Matrix- 
normen. 


Aus ||Al| < 1 folgt die Existenz [4] der Matrix 
(E — A)-! und damit die Existenz und Eindeutigkeit 
der Lösung x von (1). Wegen 


{B— Ay! AR = {H— Ar— An)! Ar 
= {{E— Ay] [B—(B-— Ar)" Anl}!Ar 
= [B—(B— Ay)" Ar "[E— Ar] Ar 


=(E—L)"!L 
existiert dann auch die Matrix (# — L)-!, es gilt 
All 
E—L)!L| = |(E — 4)" ApR| Ss -—— 7 


und nach (9) die 
Fehlerabschätzung von W. Dück: 


All 
— all s ——— _ 17): 
(m+ı z|| = 1—][A]] m+ı zm|| (17) 


(m=.0,1,2,...) 


Diese Fehlerabschätzung läßt sich, falls die ge- 
wählte Matrixnorm monoton und absolut ist, d.h. 
falls aus 


Ars |Arl + Arl= dr + Anl = 4 
die Gültigkeit von 
IAzll = II Arl || = || 1A] || = IIAll 
folgt, (vergröbern und) vereinfachen zur 
Fehlerabschätzung von L. CoLLarz: 


A 
Im +ı — ®|| sl Ian+ı mil (18). 


u a er) 


Unter der Voraussetzung ||A|| < 1 wurde die Kon- 
vergenz des Iterationsverfahrens bisher nicht unter- 
sucht; es gibt aber ein weiteres hinreichendes 


Konvergenzkriterium von J. WEISSINGER [8]: 
(19) 

Ist die gewählte Matrixnorm monoton und absolut 
und ||A|| < 1, so konvergiert die Iteration in Einzel- 


schritten (2) gegen die dann eindeutig existierende 
Lösung von (1). 


Für die nicht absolute Spektralnorm nennen wir 
ein Beispiel mit ||A||spektraı < 1, bei dem die Itera- 
tion in Einzelschritten (2) divergiert: 

1 


B: 
||A||spektraı = + ve. od)=1. 


5. Vergleich der Konvergenzkriterien und 
Fehlerabschätzungen 


Die Fehlerabschätzung (16) und das zugehörige 
Konvergenzkriterium (15) sind beide für beliebige 
Matricnormen gültig, das zur Fehlerabschätzung von 
W. Dück (17) oder L. Corarz (18) gehörige Konver- 
genzkriterium von J. Weıssınger (19) jedoch nur 
für monotone, absolute Matrixnormen. 


An Stelle der Voraussetzung yall <1l von (17), 
(18), (19) tritt bei (15), (16) die Voraussetzung 
C <1, die auch für ||A|| = 1 erfüllt sein kann. 


Die Bedingung € < 1 läßt zu, daß in der Matrix 
E — A; — Ar Elemente von 4, betragsmäßig über- 
wiegen, falls die Elemente von Ar hinreichend kleinen 
Betrages sind; insbesondere zeigt die Fehlerabschät- 
zung (16) für beliebige A, und Ar =0 die exakte 
Lösung nach einem Iterationsschritt an. 

Für ||4|| <1 kann, wie Beispiele zeigen, die 
Fehlerabschätzung (16) besser oder schlechter sein 
als die Fehlerabschätzungen von L. CoLLATZ und 
W. Dück. 


Die für die Matrixnorm der maximalen Zeilen- 
betragssumme (mit Gewichten) gültige Fehlerab- 
schätzung von H. SAssENFELD [6] (und deren Verall- 
gemeinerung von J. SCHRÖDER [7]) ist bezüglich dieser 
speziellen Norm besser oder genau so gut wie die 
Fehlerabschätzungen (16), (17), (18). 
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BUCHBESPRECHUNGEN 


W. Lietzmann, Aus meinen Lebenserinnerun- 
gen. 114$S. Göttingen 1959. Vandenhoeck & Rup- 
recht. Preis brosch. DM 6,80. 


h Die Entwicklung, die der mathematische Unter- 
richt an den höheren Schulen in Deutschland in der 
ersten Hälfte unseres Jahrhunderts genommen hat 
und die ihren entscheidenden Anstoß von FrLix KLEIN 
erhalten hat, ist zu einem bedeutenden Teil dem Wir- 
ken WALTER LIETZMANNs zu verdanken. Seine nach 
seinem Tode (1959) erschienenen, von K. FrLapr 
herausgegebenen Lebenserinnerungen lassen bei aller 
Bescheidenheit des Verfassers die außerordentliche 
Bedeutung seiner umfassenden Tätigkeit sichtbar 
werden, die die normale berufliche Arbeit eines Schul- 
mannes bei weitem übersteigt. Seine Zeit und die 
vielen hervorragenden Persönlichkeiten, mit denen er 
in Berührung gekommen ist, werden vor dem geistigen 
Auge desLesers noch einmal lebendig, und demjenigen, 
dem die heutige Problematik des mathematischen 
Unterrichts Sorgen bereitet, wird vollends klar, daß 
diese Problematik in vieler Hinsicht gar nicht neu 
ist, sondern schon von den vorangegangenen Genera- 
tionen mit großem Ernst erkannt und ihre Lösung 
auch schon von ihnen mit viel Klugheit versucht 
worden ist. 


Dresden H. HEınRıcH 


L. M. Milne-Thomson, Plane Elastic Systems. 
VIII + 221 S. m. 76 Abb. Berlin/Göttingen/Heidel- 
berg 1960. Springer-Verlag. Preis brosch. DM 49,80. 

Die mathematischen Hilfsmittel werden im Kapi- 
tel III bereitgestellt. Der Integralsatz von STOKES 
wird für die Ebene spezialisiert und für komplexe 
(nicht notwendig holomorphe) Funktionen formuliert. 
Es schließen sich an: CAaucHyvscher Integralsatz, 
CaucHvsche Integralformel, CAvcHysche Integrale 
und die Formeln von PLEMELJ. Weiterhin werden 
die HıLBertschen Randwertprobleme für Bögen for- 
muliert und die dazugehörigen stückweise holomor- 
phen Funktionen in Gestalt von Integralen hergeleitet, 
die PLEMELJ-Funktionen enthalten. 

Die allgemeine Theorie wird in den KapitelnI und II 
geboten, wo die Rolle der Aıkyschen Spannungs- 
funktion (im isotropen und anisotropen Fall) und der 
komplexen Spannungsfunktion (im isotropen Fall) 
dargelegt wird. Die hierbei benötigten Grundlagen 
der Elastizitätstheorie für ebene Systeme werden aus 
denen des räumlichen Falles spezialisiert. Diese 
wiederum werden kurz hergeleitet. Sollte hierbei eini- 
gen Lesern die benutzte Symbolik beim Gebrauch des 
Tensorkalküls (‚,‚skalare‘‘ und ‚‚doppelt-skalare‘‘ Mul- 
tiplikation von Tensoren [Bezeichnung durch 
bzw. +], Dyaden- und \V-Rechnung) ungewohnt 
sein, so ist dies nicht tragisch, weil hierauf später 
nicht zurückgegriffen wird. 

Die Kapitel IV und V beschäftigen sich mit speziel- 
len Problemen beim Vorliegen von geradlinig ge- 
schlitzter Ebene, Halbebene, Kreisscheibe, Kreisloch, 
Kreisbogenschlitzen, Kreisring und Halbkreisring. 
Zur Lösung werden die Hinserrschen Randwert- 
probleme aus Kapitel III benutzt. 

Die Zurückführung des Falles allgemeinerer Ränder 
durch konforme Abbildung auf Kreisprobleme zeigt 
Kapitel VI. 

Schließlich wird in Kapitel VII der Einfluß der 
Anisotropie diskutiert und gezeigt, daß hier keine 
grundsätzlichen Schwierigkeiten hinzukommen, son- 
dern lediglich eine Vergrößerung der Rechenarbeit 
auftritt. 

Die zu den Kapiteln hinzugefügten Aufgaben (ins- 
gesamt 150) sind nicht nur Übungsmaterial; sie be- 
leuchten den behandelten Stoff in einer vom Text ver- 
schiedenen Weise. 


Dresden H. WENZEL 


A. Duschek, Vorlesungen über höhere Mathe- 
matik, Zweiter Band: Integration und Differentia- 
tion der Funktionen von mehreren Veränderlichen. 
Lineare Algebra. Tensorfelder. Differentialgeometrie. 
Zweite, neu bearbeitete Aufl. VIII --401S. m. 


136 Abb. Wien 1958. Springer-Verlag. Preis geb. 


DM 45,—. 


Aus der Besprechung der zweiten Auflage des ersten 
Bandes (ZAMM 37 (1957), S. 159) ist bekannt, daß das 
hervorragende Lehrbuch der Höheren Mathematik, 
das 1951 seine erste Auflage erlebte und in dieser Zeit- 
schrift Bd. 31 (1951), 8. 58—59, ausführlich be- 
sprochen worden ist, eine andere Stoffeinteilung er- 
fahren hat. Der vorliegende zweite Band enthält jetzt 
die folgenden Kapitel: 

I. Grundbegriffe; Differentistion der Funktionen 
von mehreren Veränderlichen (Ergänzungen aus der 


Lehre von den Punktmengen; Funktionen mehrerer 


Veränderlicher; Grenzwert und Stetigkeit; Diffe- 
rentiation; homogene Funktionen; TAyLorsche For- 
mel; Doppelfolgen und Doppelreihen; Koordinaten- 
transformation, Punkttransformation und Abbildung 
zweier Ebenen oder Räume; ebene Kurven; Extrema 
von Funktionen mehrerer Variabler; Grundbegriffe 
der Vektorrechnung). 

II. Die Integration der Funktionen von mehreren 
Veränderlichen (Integrale als Funktionen eines Para- 
meters; Kurvenintegrale und lineare Differential- 
formen; Bereichsintegrale, mehrfache Integrale in 
Geometrie und Mechanik). 

III. Lineare Algebra (Determinanten und Matrizen; 
lineare Gleichungen; lineare Transformationen, Vek- 
toren und Tensoren; Tensoren zweiter Stufe). 

IV. Tensoranalysis und Differentialgeometrie (Be- 
griff des Tensorfeldes und Differentiation der Feld- 
größen; Integration der Feldgrößen; Raumkurven; 
Grundzüge der Flächentheorie). 

Anhang: Lösungen der Aufgaben. 

Die Vorzüge dieses Werkes sind schon oft hervor- 
gehoben worden: Es schenkt dem Leser nichts, weil es 
überall die notwendige mathematische Strenge wahrt, 
kommt ihm aber durch eine außerordentlich klare und, 
wo esangebracht ist, sehr ausführliche Darstellung ent- 
gegen. Es hält viele anschaulich faßbare Beispiele bereit 
und beachtet die Zusammenhänge mit den Anwen- 
dungen in Physik und Technik, wie es für ein Lehrbuch 
notwendig ist, das für die Studierenden unserer Tech- 
nischen Hochschulen gedacht ist. Es wäre zu wün- 
schen, daß die Zahl der wissenschaftlich arbeitenden 
Ingenieure, deren mathematische Bildung dem Inhalt 
und Niveau dieses ausgezeichneten Lehrbuches ent- 
spricht, nicht zu klein bleibt. 


Dresden H. HEINRICH 


V. Fock, Theorie von Raum; Zeit und Gravi- 
tation. XIX + 501 8. Berlin 1960. Akademie-Ver- 
lag. Preis geb. DM 42,50. 


In dem vorliegenden Buch konnte Professor Fock, 
der bekanntlich zu der Entwicklung der Relativitäts- 
theorie insbesondere in den letzten 25 Jahren Beiträge 
von grundlegender Bedeutung gegeben hat, ein nach 
zwei Richtungen zielendes Programm in besonders 
harmonischer Weise erfüllen. Erstens gibt er eine aus- 
führliche Darstellung des klassischen Inhalts des- 
jenigen Gebietes der Physik, welches man gewöhnlich 
als Relativitätstheorie bezeichnet. Darüber hinaus 
stellt er in den beiden letzten Kapiteln des Buches 
systematisch die zu einem großen Teil von ihm persön- 
lich stammenden neueren Entwicklungen in der Re- 
lativitätstheorie zusammen. Dadurch ist ein Buch 
entstanden, welches für Physiker oder Mathematiker, 
die die Relativitätstheorie kennenlernen wollen, be- 
sonders wertvoll und für den in der Relativitätstheorie 
arbeitenden Forscher unentbehrlich wird. 
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Die folgenden Zeilen mögen den Inhalt des Buches 
näher erläutern. Die Kapitel 1 und 2 behandeln die 
spezielle Relativitätstheorie in 3- und 4-dimensionaler 
Form. Das 3. Kapitel enthält eine Einführung in die 
Tensorrechnung für Rıemannsche Räume. Das 
4. Kapitel gibt dann eine Darstellung der speziellen 
Relativitätstheorie bei Verwendung von beliebigen 
Koordinatensystemen. > 

In den übrigen Kapiteln wird die Eıysreinsche 
Gravitationstheorie behandelt. Insbesondere ist in 
Kapitel 5 der klassische Inhalt dieser Theorie zu- 
sammengestellt: Die Aufstellung der Feldgleichungen 
und ihre allgemeine Diskussion, die Herleitung der 
SchwarzscHitpschen Lösung und der 3 allgemein- 
relativistischen Effekte. Die besonders tiefgehende 
Analyse des Erssreisschen Äquivalenzprinzips, die 
dieses Kapitel enthält, sei hervorgehoben. 

Kapitel 6 enthält dann eine ausführliche Darstellung 
der zu einem großen Teil vom Verfasser selbst stam- 
menden neueren Ergebnisse bezüglich des Bewegungs- 
problems in der allgemeinen Relativitätstheorie, wel- 
ches bekanntlich in der Herleitung der Bewegungs- 
gleichungen für die das Gravitationsfeld erzeugenden 
Körper aus den Feldgleichungen besteht. 

In dem letzten Kapitel wird zunächst die Diskussion 
des Bewegungsproblems weiter verfolgt, diesmal im 
Zusammenhang mit dem Problem der Gravitations- 
strahlung: Das Gravitationsfeld wird in der ersten 
Näherung für große Entfernungen von den felderzeu- 
genden Körpern bestimmt und dann die von diesen 
Körpern emittierte Gravitationsstrahlung diskutiert. 
Es folgt eine methodische Darstellung der Gedanken 
und Argumente des Verfassers, die ihn zu seiner schon 
viel diskutierten Forderung veranlaßten, die allge- 
meine Kovarianz der Gravitationstheorie mittels einer 
mit den Feldgleichungen verträglichen Zusatzbedin- 
gung — der Harmonizitätsbedingung — so einzu- 
schränken, daß man auch in dieser Theorie zu physi- 
kalisch bevorzugten Koordinatensystemen gelangt. 
Die hier gegebene ausführliche Zusammenstellung 
dieser Argumente sowie der vom Verfasser bisher er- 
reichten diesbezüglichen Ergebnisse werden den in 
diesem Gebiet arbeitenden Forschern wertvolle An- 
regungen geben. 

Es ist sehr zu begrüßen, daß dieses schon im Jahre 
1955 in russischer Sprache erschienene Buch jetzt 
auch in deutscher Sprache vorliegt. 


Berlin A. PAPAPETROU 


S. Krug — P. Stein, Einflußfelder orthogonal 
anisotroper Platten. Influence surfaces of ortho- 
gonal anisotropic plates. In deutscher und englischer 


Sprache. Englische Übersetzung durch H. Juhl. 
VIII --32 8. m. 31 Abb. u. 193 Tafeln. Berlin/ 
Göttingen/Heidelberg 1961. Springer-Verlag. Preis 


geb. DM 67,50. 


In Erweiterung des bekannten Einflußfelder-Bandes 
von PUCHER für die Schnittkräfte isotroper Platten 
werden in vorliegendem Tafelwerk die Einflußfelder 
für die Feldmitten- und Randmittenmomente der- 
jenigen orthogonal-anisotropen Platten angegeben, 
die der Husurschen Differentialgleichung genügen. 
Die theoretischen Grundlagen für die Übertragung der 
Singularitätenmethode auf orthogonal-anisotrope 
Platten und der zweckmäßige Lösungsaufbau wurden 
von den Verfassern in ihren Dissertationen erarbeitet. 
Im Sinne eines handlichen Tafelwerkes für den prak- 
tischen Gebrauch wurden in den vorliegenden Band 
von den Verfassern nur die für die Auswertung der 
Einflußfelder erforderlichen Angaben aufgenommen; 
die Auswertung wurde an Beispielen ausführlich de- 


monstriert. Die Vielzahl der aufzunehmenden Seiten- 


und Steifigkeitsverhältnisse konnte durch zwei Maß- 
nahmen auf ein erträgliches Maß beschränkt werden. 
Analog zu den Einflußflächenangaben von HoMBERG 
für die unendlich langen, an zwei Rändern frei dreh- 


bar gelagerten orthotropen Platten wurden auch hier 
die Platten. und Lastkoordinaten ‚transformiert, so 
daß zur Beschreibung der geometrischen und elasti- 


schen Eigenschaften einer Platte zwei Parameter — das 


„modifizierte‘‘ Seitenverhältnis 


I, By 


und der ‚‚effektive‘‘ Drillsteifigkeitsfaktor x — aus- 
reichten. Die Einflußfelder wurden für fünf Werte 8 
drei Werte x und sechs Randbedingungskombinatio- 
nen in 193 Tafeln angegeben. Für Zwischenwerte e 
wurden die Interpolationsmöglichkeiten erläutert, für 
abweichende <-Werte konnten Inter- bzw. Extra- 
polationsformeln angeschrieben werden. Auf die An- 
wendungsmöglichkeit dieser Einflußfelder für spezielle 
durchlaufende Systeme wurde hingewiesen. 

Neben der Möglichkeit einer einfachen und raschen 
Untersuchung der genannten (kontinuierlichen) ortho- 
tropen Platten gestattet das vorgelegte Tafelwerk 
auch Trägerroste im Sinne von Guyox und MAssoNEr 
mit guter Näherung zu berechnen. 

Es unterliegt keinem Zweifel, daß das sehr gut aus- 
gestattete Tafelwerk bald einen festen Platz in der 
Konstruktionspraxis einnehmen wird. 


Dresden G. BÜRGERMEISTER 


I. I. Gol’dman and V.D.Krivchenkov, Problems 
in Quantum Mechanics. 2758. m. 33 Abb. Ox- 
ford/London/New York/Paris 1961. Pergamon Press. 
Preis geb. 50. 


Das Buch enthält eine Aufgabensammlung mit 
Lösungen zur Quantentheorie, die für Studenten des 
4. Studienjahres gedacht ist. Die einzelnen Beispiele 
sind der praktischen Erfahrung erwachsen, die beim 
Unterricht dieses Faches an der Moskauer Lomo- 
nossow-Universität gewonnen wurden. Ihrem Thema 
nach beschränken sich die Aufgaben auf die nicht- 
relativistische Quantentheorie. Besonders sorgfältig 
sind Fragen herausgearbeitet, bei denen Drehimpuls 
und Spin eine Rolle spielen. Auch Aufgaben, deren 
Lösung das Rechnen mit Matrizen erfordert und die 
erfahrungsgemäß den Studenten besondere Schwierig- 
keiten machen, sind in die Sammlung aufgenommen 
worden. Das Buch enthält zwei Anhänge, einen über 
die WKB-Methode zur näherungsweisen Berechnung 
der Wellenfunktion, und eine kurze, aber sehr klare 
und durchsichtige Beschreibung des Isotopenspins 
von Nukleonen. 


Dresden W. MAcKE 


V. N. Faddeyeva — N. M. Terent’ev, Tables of Va- 
lues of the Function 


— z? r 2 i P t? 
w(2) = e Di = e dt 
In. 


ö 
for Complex Argument. 280 8. Oxford/London/ 


New York/Paris 1961. Pergamon Press. Preis geb. £ 5. 


Die vorliegende Tafel (eine Übersetzung aus dem 
Russischen) enthält sechsstellige Werte der obigen 
Funktion w(z) =u+iv 2e=x-+iy), die in vielen 
Gebieten der Theoretischen Physik (Ausbreitung elek- 
tromagnetischer Wellen, Quantenmechanik, Astro- 
physik u.a.) vorkommt, und zwar in Tafel I für 
© = 0(0,02)3 und y = 0(0,02)3 und in Tafel II für 
© = 3(0,1)5 und y = 0(0,1)2,9 sowie für x = 0(0,1)5 
und y = 3(0,1)5. In sieben Spalten sind x, u, v sowie 
deren modifizierte 1. Differenzen und 2. Differenzen 
(in ®-Richtung) für jeweils ein y zusammengestellt. 
Bei quadratischer Interpolation (hierfür sind geeignete 
vereinfachte Formeln angegeben) beträgt der Inter- 
polationsfehler in Tafel I höchstens ein bis zwei Ein- 


en SEE ee 0 en 2 2 in 


a na 1 in Du Lan ce 


TE 


Eingegangene Bücher 519 


heiten der sechsten Dezimale, in Tafel II — bis auf 
gewisse Bereiche in der Nähe der reellen Achse — 
ebenfalls. 

In der Einführung werden zunächst ausführlich die 
Eigenschaften der Funktion w(z) behandelt: Integral 
darstellungen, Potenzreihenentwicklung, Zusammen- 
hang mit Kettenbrüchen und asymptotisches Ver- 
halten. Da w(iy) reell ist, nimmt w(2) konjugiert 
komplexe Werte in symmetrisch zur imaginären Achse 
gelegenen Punkten an; ferner giltw(—z) =2e” "_ (2). 


Dadurch konnte die Tafel auf den Bereich > 0, 
y Z0 beschränkt werden. Dann folgen Hinweise auf 
die Berechnung der Tafel, wobei auch die durchge- 
führten Kontrollen erwähnt werden. Die Benutzung 
der Tafel wird anschließend an einigen Beispielen er- 
läutert, und zum Abschluß wird, das Berechnen von 
Werten außerhalb des tabulierten Bereichs bespro- 
chen (ebenfalls mit Beispielen). 


Dresden A. SCHUBERT 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


P. Moon and D. E. Spencer, Field Theory Hand- 
book. Including Coordinate Systems Differential 
Equations and their Solutions. VIII + 2368. Ber- 
lin/Göttingen/Heidelberg 1961. Springer-Verlag. Preis 
geb. DM 69,—. 


StructureofLanguageandits Mathematical 
Aspects. (Proceedings of Symposia Applied Mathe- 
matic, Volume XII.) VI+279S. Providence, Rh. I., 
1961. American Mathematical Society. Preis geb. 
$ 7.80. 


Selected Translations in Mathematical Sta- 
tisties and Probability, Vol. 1. V + 3068. 
Providence, Rh. I., 1961. Printed for Institute of 
Mathematical Statistiss by American Mathematical 
Society. Preis geb. $ 4.80. 


J. Aezel, Vorlesungen über Funktionalglei- 
chungen und ihre Anwendungen. 3318. m. 
14 Abb. Berlin 1961. VEB Deutscher Verlag der 
Wissenschaften. Preis geb DM 37,50. 


Zweiter Ungarischer Mathematischer Kon- 
greß, Budapest, 24.—31. August 1960. Band I 
u.1I. Budapest 1961. Akademiai Kaidö. 


S. Flügge, Lehrbuch der Theoretischen Phy- 
sik. Band I: Einführung. Elementare Mechanik 
und Kontinuumsphysik. 2568. m. 47 Abb. Berlin/ 
Göttingen/Heidelberg 1961. Springer-Verlag. Preis 
geb. DM 25,80. 


6. B. Thomas, Calculus. Sec. Ed. XIII + 850 8. 
m. 348 Abb. u. 4 Tafeln. London 1961. Addison- 
Wesley Publishing Comp. Preis geb. $ 8.75. 


S. Vajda, Mathematical Programming. IX + 
3108. London 1961. Addison-Wesley Publishing 
Company. Preis geb. $ 9.00. 


E. Pestel, E. Kollmann, Grundlagen der Rege- 
lungstechnik. Ein Lehrbuch für Studierende und 
Ingenieure. (Regelungstechnik in Einzeldarstellungen, 
Band 1.) VII-+ 3228. m. 396 Abb., 20 Tab. u. 
148 Übungsaufgaben. Braunschweig 1961. Friedr. 
Vieweg & Sohn. Preis geb. DM 33,—. 


Fr. A.Willers, Elementar-Mathematik. 10. Auf- 
lage. VIII -+ 267 S. m. 172 Abb. Dresden und Leip- 
zig 1961. Verlag von Theodor Steinkopff. Preis geb. 
DM 16,—. 


A. Haerter, Theoretische und experimentelle 
Untersuchungen über die Lüftungsanlagen 
von Straßentunneln. (Mitteilungen aus dem In- 
stitut für Aerodynamik an der Eidgenössischen Tech- 
nischen Hochschule in Zürich, Nr. 29.) 1018. Zürich 
1961. Verlag Leemann. 


N. Athanassiadis, Potential Flow through Spi- 
ral Casings. (Mitteilungen aus dem Institut für 
Aerodynamik an der Eidgenössischen Technischen 
Hochschule in Zürich. Nr. 30.) 1758. Zürich 1961. 
Verlag Leemann. 


B.M. Brown, The Mathematical Theory of 
Linear Systems. XV + 2678. London 1961. 
Chapman & Hall Ltd. Preis geb. 50 s. net. 


Jordan-Eggert-Kneisel, Handbuch der Vermes- 
sungskunde. Zehnte, völlig neu bearbeitete und 
neu gegliederte Ausgabe. BandI: Mathematische 
Grundlagen, Ausgleichsrechnung und Rechenhilfs- 
mittel. XVI-+ 8088. m. 208 Abb. Stuttgart 1961. 
J. B. Metzlersche Verlagsbuchhandlung. Preis geb. 
DM 160,—. 


V.L. Zaguskin, Handbook of Numerical Me- 
thods for the Solution of Algebraicand Trans- 
cendental Equations. Translated from the Rus- 
sian. XIX +195S. Oxford/London/New York/ 
Paris 1961. Pergamon Press. Preis geb. 40 s. 


N. A. Archangelski — B. I. Saizew, Automatische 
Ziffernrechenmaschinen. 1308. m. 23 Abb. 
Berlin 1960. VEB Deutscher Verlag der Wissen- 
schaften. Preis brosch. DM 5,80. 


$. Hartman and J. Mikusinski, The Theory of 
Lebesgue Measure and Integration. Translated 
from Polish. 176 S. Oxford/London/New York/ Paris 
1961. Pergamon Press. Preis geb. 30 s. 


H. Benz, Über eine Bewertungstheorie der 
Algebren undihre Bedeutungfür die Arithme- 
tik. (Heft9 der Schriftenreihe der Institute für 
Mathematik.) 149 S. m. 5 Abb. Berlin 1961. Akade- 
mie-Verlag. Preis brosch. DM 29,50. 


H.F. Chilmi, Qualitative Methoden beim 
n-Körperproblem der Himmelsmechanik. 
(Heft 10 der Schriftenreihe der Institute für Mathe- 
matik.) 1168. m. 4Abb. u. 1 Tab. Berlin 1961. 
Akademie-Verlag. Preis brosch. DM 23,—. 


W. Klein, Grundlagen der Theorie elektri- 
scher Schaltungen. XII+ 2598. m. 120 Abb., 
3 Kunsttafeln, 1 Übersichtstafel u. 2 Tab. Berlin 
1961. Akademie-Verlag. Preis geb. DM 39, —. 


0.-H. Keller, Zur Berechnung desGalois’schen 
Körpers und der Galois’schen Gruppe einer 
Gleichungin erträglich vielen Schritten. (Be- 
richte über die Verhandlungen der Sächsischen Akade- 
mie der Wissenschaften zu Leipzig, Mathematisch- 
naturwissenschaftliche Klasse, Band 104, Heft 5.) 
218. Berlin 1961. Akademie-Verlag. Preis brosch. 
DM 2,—. 
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G.Rackwitz, PraktischeMathematik. (Taschen- 
bücher für Holztechnik und Holzwirtschaft, Bd.T.) 
XII+ 1478. m. 55Abb. Prien/Chiemsee 1961. 
Helmut Bücking Verlag. Preis brosch. DM 14,80. 


H. Hasse / W. Klobe, Aufgabensammlung zur 
Höheren Algebra. 3. Aufl. (Sammlung Göschen, 
Band 1082.) Berlin 1961. Walter de Gruyter & Co. 
Preis brosch. DM 3,60. 


F. Sehytil, Wirbelschichttechnik. VII + 
120 S. m. 87 Abb. Berlin/Göttingen/Heidelberg 1961. 
Springer-Verlag. Preis brosch. DM 16,50. 


H. Hermes, Aufzählbarkeit, Entscheidbar- 
keit, Berechenbarkeit. Einführung in die Theorie 
der rekursiven Funktionen. X -+ 2468. m. 3 Abb. 
Berlin/Göttingen/Heidelberg 1961. Springer-Verlag. 
Preis geb. DM 49,80. 


D. Morgenstern und I. Szabö, Vorlesungen über 
Theoretische Mechanik. XII + 374 S. m. 
112 Abb. Berlin/Göttingen/Heidelberg 1961. Sprin- 
ger-Verlag. Preis geb. DM 69,—. 


S. Glasstone und Mc. Bilung, r \ 
theorie. Eine Einführung. X (1 S. m. 
Wien 1961. Springer-Verlag. Preis geb. DM £ 
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* @. Lefort, Algebre et analyse exercises. XI+ 
Preis 


5178. m. 29 Abb. Paris 1961. Dunod. 
44 NE. „40 I: ei 
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H. Wolt, Ausgleichsrechnung nach der 
Methode der kleinsten Quadrate. 3. Lieferung. 


Hamburg 1961. Hanseatische Verlagsanstalt GmbH. 
ee 


Preis brosch. DM 6,80. 

V.M. Blanco and $. W. MeCuskey, Basic Physies 
of the Solar System. XI + 3078. - London 1961. 
Addison-Wesley Publishing Comp. Preis geb. $ 8.50. 


W.T. Martin and E. Reissner, Elementary Dif- 
ferential Equations. Sec. Ed. XIII + 3318. 
London 1961. Addison-Wesley Publishing Comp. 
Preis geb. $ 6.75. 


W. Furrer, Raum- und Bauakustik, Lärm- 
abwehr. 2588. m. 195 Abb. Basel und Stuttgart 
1961. Birkhäuser Verlag. Preis geb. sFr. 38.50. 


NACHRICHTEN 


Internationales Kolloquium über Anwendungen 
der Mathematik in den Ingenieurwissenschaften 


In der Zeit vom 8.—10. Juni 1961 führte das Insti- 
tut für Mathematik der Hochschule für Architektur 
und Bauwesen Weimar unter Leitung von Herrn 
Prof. Dr. habil. H. MATzke ein Internationales Kollo- 
quium über Anwendungen der Mathematik in den 
Ingenieurwissenschaften mit dem Rahmenthema 
„Anwendungen elektronischer Rechenanlagen in der 
Bau- und Baustoffindustrie‘‘ durch. Die Veranstalter 
hatten es sich zur Aufgabe gemacht, über dieses aktu- 
elle Thema erstmalig in der DDR einen breiten Erfah- 
rungsaustausch anzuregen. Teilnehmer waren etwa 
150 Mathematiker und Ingenieure aus beiden Teilen 
Deutschlands, aus der UdSSR, CSSR, aus Ungarn, 
Polen und Indien, die sich mit Problemen der Bau- 
mechanik und der Ökonomie des Bauwesens be- 
schäftigen. 


Folgende Hauptvorträge wurden gehalten: 

S. Fark, Braunschweig, Knickbiegung von Durch- 
laufträgen; A. June, Jena, Ingenieurtechnische 
Probleme im Rechenzentrum Jena; H. KRETZSCH- 
MAR, Dresden, Einsatz elektronischer Rechenanlagen 
im Stahlbau; M. Rözsa, Budapest, Berechnung von 
Tragwerken mit Hilfe von Digitalautomaten; W. ZER- 
nA, Hannover, Betrachtungen und Erfahrungen über 
den Einsatz elektronischer Rechenanlagen im Bau- 
ingenieurwesen; R. ZURMÜHL, Darmstadt, Anwendung 
der Matrizenrechnung in der Statik. 

Weiterhin informierten 18 Kurzvorträge u. a. über 
den Einsatz von Ziffernrechnern und Analogrechnern 
im Bauwesen sowie über Anwendungen der linearen 
Programmierung in der Bauökonomie und Gebiets- 
planung. Die Vorträge sollen in einem Sonderheft der 
Wissenschaftlichen Hochschulzeitschrift mit sämt- 
lichen Vorträgen herausgegeben werden. 

o. H. MATZKE 


BERICHTIGUNGEN 


Zu H. Buchholz, Die Erwärmung in metallischen 
Kreislochscheiben .. ., ZAMM 41 (1961), S. 229—238: 
In der obigen Arbeit sind folgende Änderungen vor- 
zunehmen: 
1. Gleichung (2) auf $. 231 lautet: 
d2d(vs; pP) 1 ddw;P), - v2 
”. avi vo dvo Pe; Pl =) 
Op |] - 
a =. 80 > a): 
Da nl=se ..@) 
2. In Gleichung (7) auf 8.232, 2. Zeile, lautet im 
Argument von Jy der erste Faktor i b anstelle von b. 
3. Auf 8.233 ist in der Formelzeile zwischen den 
Gleichungen (3) und (3a) o durch 02 zu ersetzen. 


4. Auf S. 234 ist in den beiden Textzeilen unter 
Gleichung (6a) y durch » zu ersetzen. 

5. Auf 8. 236 ist in der zweiten Zeile der Gleichung 
(14) die letzte Klammer zu streichen. 


6. Auf S. 237 steht in Gleichung (18b) h?,, anstelle 
von Ay,r- 


Zu B. Tanimoto, The Reciprocal Theorem for 
Applied Mechanics, ZAMM 41 (1961), S. 185—197: 


Der Verfasser teilt mit, daß auf S. 197 in dem Aus- 
druck für w’ auf der vierten Zeile das letzte Plus- 
Zeichen zu streichen ist. 


Herausgeber und Haupbschriftleiter: Prof. Dr.-Ing. habil. H. Heinrich, Dresden A 27, Friedrich-Hegel-Str. 31. Verlag: Akademie-Verlag 
GmbH, Berlin W 8, Leipziger Straße 3-4; Fernsprecher: 220441, Telex-Nr. 011773, Postscheckkonto: Berlin 35021. Bestellnummer dieses 
Heftes: 1009/41/12. Die Zeitschrift für Angewandte Mathematik und Mechanik erscheint monatlich. Bezugspreis: vierteljährlich DM 15,—. 
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WISSENSCHAFTLICHE GRUNDLAGEN 
DER MODERNEN TECHNIK 


Herausgegeben von Prof. Dr. Hans Frühauf, Prof. Dr. Robert Rompe, 
Prof. Dr. Kurt Schröder, Prof. Dr. Erich Thilo, Prof. Dr. Hans Gummel, 
Prof. Dr. Hermann Klare, Prof. Dr. Eberhard Leibnitz 


Reihe A: Tagungen 
Band 1 


ELEKTROTECHNIK UND MASCHINENBAU 


Vorträge, gehalten auf der Tagung der Forschungsgemeinschaft 
am 8. und 9. 12. 1960 in Berlin-Adlershof 


1961. 168 Seiten — 67 Abbildungen — 24 Tabellen — gr. 8° — Halbleinen DM 9,80 


Am 8. und 9. Dezember 1960 veranstaltete die Forschungsgemeinschaft der naturwissenschaft- 
lichen, technischen und medizinischen Institute der Deutschen Akademie der Wissenschaften zu 
Berlin eine Konferenz ‚‚Über wissenschaftliche Grundlagen der modernen Technik — Elektro- 
technik und Maschinenbau“. 

Die Ergebnisse der Konferenz liegen mit dieser Arbeit vor. Sie dienten der Entwicklung und 
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Festigung der wissenschaftlichen Grundlagen für die Produktion. Namhafte Wissenschaftler und 
Praktiker berichten über Entwicklungstendenzen auf dem Gebiet der Elektrotechnik und im 
Maschinenbau, über Probleme moderner physikalischer Analysenmethoden, über Gebiete neuer 
metallischer Werkstoffe für Elektrotechnik und Maschinenbau und über Kunststoffe mit be- 
sonderen physikalischen Eigenschaften. 


. Die Vorträge geben in einer leicht verständlichen wissenschaftlichen Darstellung einen ausge- 


zeichneten Überblick über den Stand und die Entwicklungsrichtung des jeweils behandelten 
Gebietes. Sie zeigen die Bedeutung, die die Anwendung.der neuesten Forschungsergebnisse für 
die industrielle Fertigung besitzt. Deshalb stellt sich die Forschungsgemeinschaft mit dieser 
Schriftenreihe die Aufgabe, diese Ergebnisse aus dem Bereich der naturwissenschaftlich-tech- 
nischen Institute der Deutschen Akademie der Wissenschaften zu Berlin schnell unserer In- 
dustrie zur Verfügung zu stellen. Damit soll erreicht werden, daß unsere industriellen Produkte 


den Höchststand der modernen Technik aufweisen. 


In Vorbereitung 


Band 2 
KONTAKTWERKSTOFFE IN DER ELEKTROTECHNIK 


Bestellungen durch eine Buchhandlung erbeten 


AK CA DUE-M I EN. : BER LIN 


Hera RUE den Ha Professoren Dr 
Mol a Dr. W. DEKEYSER, ( u 
Hull — Dr. Dr. h. c. P. GÖRLICH, Jena_Dr. E | 
Dr. P. T. LANDSBERG, Cardiff — Dr. L. N “ e—: 
- SEEGER, Stuttgart — Dr. O. STASIW, Berlin — D dr. M. 
MANN, Karlsruhe — Dr. 6. SZIGETI, Budapest unc 


in Taf an 


Schriftleitung: Prof Dr. PAUL GÖRLICH, Jena ”E = 3 | 
Redaktionskollegium: : Dr. 8. OBERLÄNDER, Drn.Gi 


Jen 10 Shen mi Alumni Format 624 em, um P 
Erscheinungsweiset monatlich E € : 


In dieser neuen Zeitschrift Se Arbeiten der Grundlagen- 
" forschung aus dem Gesamtgebiet der sie ver- 
öffentlicht. - 


Insbesondere wird über ip Atommistik der Meaikkiepeniehle 
schaften, zum Beispiel über Bau und BURIDE von Kristallen 
Die monatlich erscheinenden Hefte bestehen aus re Pe 
Zusammenfassende Berichte über aktuelle Teilgebiete der Festkörperphysik 3 
Originalarbeiten, die nicht schon an anderer Stelle veröffentlicht wurden und neue w 


Erkenntnisse enthalten 
Kurze Originalmitteilungen 


Vorabdruck des Titels und der kurzen Zusammenfassung von Örizinalarheiige dieser 7 
(in englischer Phersetzung) der „‚Fizika Tverdogo Tela“ 


Die Zeitschrift physica status solidi soll als internationales 
Organ der Festkörperphysik eine möglichst rasche Veröffent- 


dern, die bestehenden 
gänzen. Die Arbeiten erscheinen in deutscher, e 


lichung der neuen Erkenntnisse und Forschungsergebnisse 
garantieren, und durch Konzentration wichtiger Arbeiten 
dieses Fachgebietes, insbesondere aus den europäischen Län- 


Bestellungen durch eine Buchhandlung erbeten 


AKA DE MI E- 


VERLIESS 


französischer oder russischer Sprache, Eine kurze Zı N 
fassung der Arbeiten wird auch in einer en Sea 
der des 1 a Di 


